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1.0 INTRODUCTION 

This task is the fourteenth in the series accomplished within the 
scope of the Work Statement, Attachment A, Item 5 of Contract 
NAS5-23185 . 


2.0 SCOPE 

This task was issued to perform the following analysis and asso- 
ciated effort: 

Perform suitable analysis, assuming Lambertian surface, to deter- 
mine required sensor offset in a lunar tracking rocket flight. 

The offset is defined as the angular difference, in arc minutes, 
between the geometrical center of the Moon and the sensor null. 
The following shall be furnished: 

• A computer program which will provide the offset, in arc 
minutes, as a function of the following lunar variables. 

a) Fraction illuminated 

b) Apparent size of Moon 

• A technical note documenting the analysis and describing 
the computer program. 


3.0 FINANCIAL SUMMARY 

The following is the final manhour and cost totals for this 
assignment . 




Labor Dollars 

MODC 

Total Costs 


Manhours 

(Unburdened) 

(Unburdened) 

W/O Fee 

al Budget 

295 

2,120 

459 

6,482 

ual Total 

361 

2,654 

433 

8,549 


4.0 CONCLUSIONS 

The analysis and results are summarized in the Final Technical 
Report F74-13 dated 31 December 1974 attached as Appendix A. 


The overrun in the manhours and total cost were primarily due to 
more than estimated time to complete the final report and the tests 
requested by GSFC that were not included in the original scope of 



the program. The tests consisted of moon intensity measurements 
made using the lunar sensor to observe the moon. This data was 
used to calibrate the sensor for the actual launch. 
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Section I 
INTRODUCTION 


This final report is submitted to Goddard Space Flight Center 
in fulfillment of a task intended to allow the center of 
brightness of the lunar crescent to be computed. This capa- 
bility was required to support sounding rocket launches from 
Kauai, Hawaii (3 November 1974) and White Sands Missile Range 
(28 December 1974) . 

This paper discusses, briefly, the operational characteristics 
of the lunar sensor which was used to point the sounding 
rocket, develops the associated mathematical model of the 
system, describes the computer programs which were written 
to implement the model, and presents data pertinent to the 
two launches. 

The work was performed under NASA Contract NAS5-23185, Task 
Assignment 014. 
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Section II 


DESCRIPTION OF THE PROBLEM 


An angular error sensor has been built by Ball Brothers Research 
Corporation for the purpose of pointing a sounding rocket pay- 
load at the moon. The sensor has an essentially linear trans- 
fer characteristic in the vicinity of null and angular position- 
errors of the vehicle are manifested as proportional error sig- 
nals which can be used by the rocket control system to reduce 
these errors to zero. 

The application for which the sensor was designed necessitated 
pointing the payload at the center of the moon; that is, the 
center of the disk which would be visible from the 'ea rth at 
full moon. 

The sensor is an electro-optical device of the "energy-balance" 
variety. Such sensors split their fields of view into two parts 
and generate error signals which are proportional to the dif- 
ference of energies gathered by each half. Thus, such a sensor 
produces a zero error signal when it is pointed at the center 
of intensity of the source. 

For the case at hand the source is, of course, the moon and for 
conditions other than full moon, the center of intensity will 
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not coincide with the selenographi-c center. It Is necessary 
to know just where the sensor will be pointing when it is nulled 

1 

so appropriate biases can be applied to correct the situation. 

The problem, then, is to determine where the center of intensity 
of the visible crescent of the moon is with respect to the 
selenographic center. In order to compute this quantity, it 
is necessary to know the realtive positions of the sun, moon, 
and sensor, the reflective characteristics of the lunar surface, 
and the optic-operational basis of the sensor. 

This report is devoted to constructing a mathematical model of 
the situation, discussing the computer program which implements 
the model, and presenting results applicable to actual flights-. 
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Section III 
MATHEMATICAL MODEL 


3.0 ASSUMPTIONS 

The simplifying assumptions upon which the mathematical model 
is based are as follows: 

1. ' The moon's albedo is not a function of wavelength. 

That is, the incident and reflected energies have 
the same spectral distribution. 

2. ’ The moon's albedo is not a function of seleno- 

graphic latitude and longitude. 

3. The sun is considered, for illumination purposes, 
to be an infinitely distant point source. Thus, 
the radiant energy is constant in the vicinity of 
the moon and,, since there is no penumbra, the 
terminator is sharply defined. 

3.1 Basic Concepts 

The lunar sensor is a two axis device and, as such, produces 
error signals about its yaw and pitch axes. To see how these 
signals are generated, consider the operation of just one of 
the channels; e.g., yaw. 
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A plane may be associated, ficticiously , with the yaw channel. 
This plane is defined by the yaw axis and optical axis of the 
sensor and it contains these axes. For want of a better term, 
we can call it the "yaw-null plane". When the sensor is pointed 
at a source of radiant energy, this plane effectively divides 
the source into two parts. The part lying on the "left" side 
of the plane is converted into one proportional voltage while 
the part lying on the "right" side of the plane is converted 
into another proportional voltage. These voltages are dif- 
ferenced to produce the final error signal which is used by 
the control system. 

Exactly the same situation applies to the orthogonal pitch axis 
and the origin of the term "energy balance sensor" becomes ob- 
vious. This concept of the sensor's operation is essential 
to understanding why the mathematical model is constructed the 
way it- is . 

3.2 Coordinate Systems 

Figure 3-1 shows the basic elements of the system; the sensor 
together with its yaw-null plane, the moon, a coordinate system 
(with axes labeled 1, 2, 3) with its origin at the center of 
the moon, and a vector, S, pointing towards the sun. 
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Figure 3-1 

Pictorial View o£ Sensor-Moon Coordinate System 


The sensor is located on the 2 -axis and the sun lies in the 
2-3 plane. These two facts serve to define the orientation 
of the coordinate system. Although it will be convenient to 
refer to the 2-3 plane as the "equator" and the 1-axis as the 
"north pole", it must be clearly understood that the true 
lunar equator and pole are something else entirely and are 
of no concern in the present development due to assumption 
2 (Section 3.0). Thus we shall use the terms equator and pole 
freely in the context just described. 

Figure 3-2 is a top view of Figure 3-1, includes no additional 
information, and is included only to clarify the situation. 

The sensor's optical axis lies in the 2-3 plane at all times 
and, as shown in Figure 3-1, the yaw-null plane intersects the 
moon somewhere between the pole and its eastern limb. This 
line of intersection, which lies on the moon's surface, is 
called M. 

Nestled within the positive 1, 2, and 3 axes is octant 1 of 
the lunar sphere; it is outlined in Figure 3-1. This octant, 
the coordinate system, and some quantitative terms are shown 
in Figure 3-3. For instance, there is the vector S again, still 
lying in the 2-3 plane, and pointing at the sun. Note that 
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Figure 3-2 


Top View o£ Sensor-Moon Coordinate System 





it is displaced an angle a from the 2-axis. As shown, a is 
positive, as are all the angles shown in Figure 3-3. 

The intersection line, M, also appeareth. It intersects the 
equator an angle 0o from the 2-axis. It is drawn so that M 
is parallel to the 1-2 plane and this is a departure from 
reality in that it results from a plane which makes a small 
non-zero angle with the 1-2 plane. This angle is never greater 
than 16 arc-minutes however (or the yaw-null plane would fail 
to intersect the moon at all) so it has been ignored. The 
resultant errors are' very small compared to those introduced 
by the assumptions discussed earlier. 

3,3 Construction of the Mathematical Model 

Our ultimate goal is to determine just where the lunar sensor 
will attain a null condition. This is equivalent to calculat- 
ing, for a given value of a (sun position) where the intersec- 
tion line, M, must be located so that the energies reaching 
the sensor from both sides of it are equal. 

To this end, we subdivide the portion of the moon visible from 
the sensor into differential areas , dA, compute the relative 
energy reaching the sensor from each of them, and then sum 
these differential energies to obtain the net effect. By com- 
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paring the energy contributed by the area to the west of M to 
the total energy, we can determine the value of 0 q which makes 
this ratio equal one-half. This, then, will be the location 
on the equator where the sensor will be nulled. 

It is now time to get down to specifics. If you are not in- 
terested in getting involved in the mathematics, this is the 
time to read elsewhere. Appendix B comes highly recommended. 

Figures 3- 4a and 3-4b are somewhat more detailed versions of 
Figure 3-3 in that they locate a typical differential area, 
dA, on the lunar surface. The manner in which dA is defined 
in each case is different and two figures are presented to un- 
clutter the artwork. 


In Figure 3-4a, dA is defined in terms of its normal vector, 

N. This vector is situated an angle <j> above the equatorial 

plane and its projection in that plane makes an angle 0 with 

1 

the 2-axis. Thus, N can be written as 


N 


s<j> 

C0C<j> 
s 6c<f> 


3.1 


^For an exposition of the perhaps unfamiliar notation used in 
this paper, see Appendix A. 
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Note that the moon is assumed to be spherical and for now, of 
radius 1 . 


Since the sensor responds to energies which are distributed 
symmetrically about the line M, we must also locate dA in a 
way that will mathematically reflect this symmetry. This is 
done in Figure 3-4b. Our differential area is now defined by 
the vector V which has its "tail" located on the 3-axis, lies 
in the sensor’s yaw-null plane, and is an angle a above the 
equatorial plane. It is written as 

V = 


c8osa 
c0 o ca 
0 


3.2 


Note that V is not, in general, a unit vector. V and N are 
related by the following equation: 


N 


0 

0 

S0o 


+ V 


3.3 


Some of the sunlight incident upon dA will ultimately be re- 
ceived by the lunar sensor. Let us now compute just what the 
intensity of this reflected sunlight will be. 
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Toward this end we let: 

B q = Power density o£ sunlight in the vicinity o£ the 
moon. (watts-cm 2 ) . 

B = Power density of the reflected sunlight at the 
point of observation. 

Y = Angle between the solar vector, S, and N. 

£ = Angle between the 2-axis (also called E) , and N. 

R = Distance from dA to the point of observation. 

f(y,£) “ a "bidirectional reflectance function" which is 

the ratio of reflected- to-incident power den- 
sities at dA. The following two attributes of 
f should be noted: 

1. It is a function of the incident angle y, 
and the reflectance angle, only and does 
not depend on the location of dA. This is 
a result of assumption 2, Section 3.0. 

2. The reflectance angle is taken to be £ rather 
than the angle between N and the line con- 
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necting dA and the observation point. This 
will produce, at worst, a 16 arc-minute error 
in the reflectance angle, the implications 
of which are so small so as to make the 
simplification more than reasonable. 

The differential power density of sunlight, dB , which emanates 
from dA and is received by the sensor is then 


dB = kB dA 

° R 2 


3.4 


where k is simply a proportionality constant that will disap- 
pear shortly when we come to compare energies reciived from 
two areas on the moon. 


The problem now is to find, for a given value of cr , the Go 

which causes line M to divide the total energy received from 

the moon into two equal parts. That is, if B(Q 0 ) is defined 

to be the result of sunlight reflected from the visible region 

of the moon to the west of M, then we seek that particular 
* 

9 0 , 0o» that will give 


B (9 o ) 

BO/2) 


1/2 


3.5 
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To accomplish this , we will have to integrate dB over a and 
0 o • 

From Figure 3- 4b, we see that 

dA = r 2 c0 o da d0 o 3.6 

where r is the radius of the moon. 

Let d be the distance from the center of the moon to the ob- 
m 

servation point. It is 'measured along the 2-axis. From the 
law of cosines we have 

R 2 = r 2 + d 2 - 2rd c£ 3.7 

m m 


Now 


c£ = N-E 


3.8 


= c<j> C0 


3.9 


so 3.7 becomes 


R 2 = r 2 + d 2 - 2rd c<f> c0 
m m 


3.10 
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Combining 3.4, 3.6, and 3.10 together with the definition 


q 


£ 



3.11 


we have 


dB = kB q : 


f C y,5)c9o 


l-2q c<j> c0 + q 5 


dad0 o 


3.12 


To actually perform the calculations, we must have y, £ , <f>, 
and 0 in terms of a and 0 O * By expanding 3.3 we get 


so 


and 


S0 


c9 osa 


C^)C0 

= 

c0 o ca 

3.13 

C<(> C0 


S0 O 


4* * 

. -1 
sm 

| c0 o sa 

| 3.14 

0 = 

-1 

tan 

I tan 0o 
ca 

| 3.15 


3-15 




The reflectance angle, £, can be obtained by combining 3.9, 
3.14, and 3.15. Finally, we get the incidence angle, y, by 
noting that 


cy = N*S 


3.16 


Now 



3.17 


so 


N*S = co c0 c4> + so s0 c<f) 


3.18 


and 


cy = c<)> c(0-o) 3.19 

which, when combined with 3.14 and 3.15 will yield y. 

We now have only to settle on what the limits of integration 
shall be. We can choose between two approaches. 
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1 . 


Integrate over the region 


- TT / 2 < a < + TT / 2 

A 

-ir/2 < 0 O 5 6o 

and construct £ so it will be zero in those 
regions not illuminated by the sun; or 

2. Integrate only over the region which is visible 
from the observation point and illuminated by 
the sun. 

The latter approach is definately preferable because it in- 
volves significantly less computation time; an important con- 
sideration since numercial integration on a digital computer 
is the only reasonable way to integrate 3.12. 

Before proceeding, we should note that because of the inherent 
symmetry of the problem, we can confine our calculations to 
the northern hemisphere and then simply double the results. 
Even the doubling can be omitted because of the ratioing 
(equation 3.5) that will eventually be done. 
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In order to integrate only over the illuminated protion of the 
moon's surface, we must be able to mathematically define the 
terminator in terms of 0o and a. Assumption 3 permits us to 
do this with a minimum of fuss. 

The terminator is defined to be the locus of points on the 
surface of the moon where S is tangent to the surface; equi- 
valently, where S and N are perpendicular. From 3.19, we see 
that this condition is met when 

c<}> c(0-a) = 0 3.20 

Now this equation is satisfied when either 


4> = ± ir/2 


3.21a 


or 


0 = o ± ir/2 3.21b 

All 3.21a states is that the terminator passes through the 
north and south poles ; a true thing but of no particular 
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interest. It is equation 3.21b that will provide a useful 
relationship . 

Before developing 3.21b, let us agree to restrict the analysis 
to values of 0 in the range 


0 < a < TT 


3.22 


This will cover the range from full moon to eclipse but will 
always produce a center of brightness in the eastern half of 
the moon. This does not cause a problem, however, since the 
basic symmetry of the problem allows us to simply negate our 
result (placing the center of brightness in the western half) 
if the original a should happen to lie in the range 


-it < 0 < 0 


3.23 


Now, with a restricted according to 3.22, we can remove the 
amgiguity of 3.21b and rewrite it as 

0 = o- tt/2 3.24 
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Combining 3.24 and 3.15 gives 


tan 


tan 6 o 


ca 


= ct-tt/2 


3.25 


Solving 3.25 for a gives the particular value of a which will 
place V on the terminator for the specified 6 0 . This special 
value of a will be called a 0 . Thus 


a 0 


cos 


-1 


■tan a tan 8 0 


3.26 


Notice that 3.26 is not defined for some combinations of o 
and 0o. This does not worry us, however, since those cases 
for which 


| tan o tan 0 O ! > 1 3.27 

are precisely those for which V is in the dark portion of the 
moon regardless of a and we will not be integrating over those 
regions anyway. 

The problem must now be separated into two basic cases depend- 
ing on the location of the sun; 
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Case I 


0 < a < ir/2 


Case II ir/2 < a < it 

The reason for this can be best understood by referring to 
Figures 3-5a and 3-5b which depict typical situations for 
Case I and Case II respectively. 


1 

TERMINATOR i 


I ' 

t TERMINATOR 






First consider 3- 5a which corresponds to Case I. The sketch 
is drawn for a a of approximately tt/4. The terminator inter- 
sects the equator at 0 o - o-tt/2 so we begin integrating there 
and let 6 0 cover the range 


o-tt/ 2 < 0 O < ir/2 

Now, over the semi-range o-tt/2 < 0o < 0, the terminator is 
visible so a will take on values 0 < a < ao where ao is cal- 
culated from 3.26. However, for positive values of 9 0 (0o in 
the eastern hemisphere) the terminator disappears onto the back 
side of the moon so a must cover the region from equator to 
limb ; i.e., 0 < a < it/ 2. 

For Case II (Figure 3-5b) the range of 0 O is constant and 
independent of a; 0 < 0o < tt/2. Again we have two subcases 
for the a interval. From Figure 3- 5b they can be seen to be: 

for 0 < 6o < o-tt/2 , a 0 < a < tt/2 

for o-tt/2 < 0 O 5 tt/2 , 0 < a < tt/2 
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3.29 

3.30 
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So there we have it. If we define the function A to mean 
"the applicable range of...", then we can write our basic 
integral as 


B(o) 



jCLiil 

l-2q c(f> c6 + q : 


da d0 o 


A (0o) A (a) 


where k is a factor of proportionality, and A (a) and A ( 0 0 ) 
are defined in the preceding paragraphs. 

We have yet to get specific about f(y,£) and come up with a 
means of calculating o. Read on. 

3.4 The "Reflectance Function", f(y,£) 

At the beginning of this study, the author (Rose) indicated 
to GSFC his intent to model the reflective properties of the 
lunar surface after a Lamertian reflector; that is 

f(y,C) = cy c? 


3.31 
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It was common "knowledge" by workers in the field that this 
was a reasonable thing to do. Nevertheless, as it was decided 
to perform some calculations using this model and compare the 
results to empirical data, thereby verifying the model's ap- 
plicability . 

The inappropriateness of the model was demonstrated instead. 

The analytical test which produced these disturbing results 
was a computation of the ratio of half moon brightness to full 
moon brightness. The Lambertian model produced (see Appendix 
C) a value of 0.32 while actual measurement' 1 ' of the ratio 
gives of value of 0.089. 

The wide discrepancy between theoretical and measured results 
initiated a search for a better reflectance function. This 
led to discussions with personnel at the Lunar and Planetary 
Laboratory, University of Arizona, in Tucson. They suggested 
that we model the lunar surface as a Lommel-Seeliger reflector; 
that is , have 


i 

C.W. Allen, ' Astrophsical Quantities', 2-nd Edition University 
of London, the Athlone Press 1964; Table of Moon's Phase Law - 
Page 146 
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f ( y, £) 


cy 

cy + c£ 


This did not help the situation, however, for when the half- 

to-full moon brightness ratio was computed, we obtained a 

value of 0.5. This was worse than the Lambertian model! It 

was then suspected that the University of Arizona information 

was being misinterpreted and that equation 3.32 gave the ratio 

of reflected to actual incident energies. If that were true, 

3.32 would have to multiplied by cos y to account for the non- 

/ 

perpendicularity of dA and S. This would give 

ffy £1 = — — 

C y + c£ 

But alas. When the half-to-full moon brightness ratio was cal- 
culated using 3.33 the situation became even worse; the value 
was 0.62. 

By now, the sounding rocket launch was becoming imminent and 
a usable reflectance function had to be settled upon soon. We 
therefore returned to the best function we had (the Lambertian 
model, equation 3.31) in hopes that even though it produced 
invalid brightness ratios, it might permit the significantly 
different problem of locating sensor null points to be solved. 


3.32 


3.33 
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To see that this might be possible, consider the £ollo\^ing. 
Suppose we sketch the lunar albedo energy falling to the west 
of the sensor yaw null plane as the plane is scanned from west 
to east across the moon. Two such plots are shown in Figure 
3-6 and they correspond to two different values of a. 



1.0 


0o 


Figure 3-6 

Sketch of B(0 O ) for Full and Half Moons 
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Let us now sketch another pair o£ simular curves only this time 
we hold o constant at tt/2 (half moon) and plot the- real B(0 O ) 
and the B(0 O ) corresponding to a Lambertian moon. 


B{0o) 



Figure 3-7 

B C 0 o ) Curves for Real and Lambertian Moons 

The point is that the sensor will point to that value of 0o 
where the B(0 O ) curve has half its end point (0 O = tt/2) value 
and that even though the two curves shown in 3-7 differ in mag- 
nitude by more than a factor of 3, their mid-points could occur 
at almost identical places. 

It was on the basis of such reasoning that the sensor offset 
was measured one night (8 h U.T., 8 September 1974) by using 
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the sensor itself and the results compared to theoretical pre- 
dictions. There was agreement to within better than 11 arc- 
seconds; much better than the goal of ±1 arc-minute for which 
we were striving. 

Thus, as a result of this very close agreement between theory 
and experiment, it was decided to use the Lambertian model in 
all subsequent calculations. It is nevertheless most annoying 
that we cannot calculate absolute brightnesses and we plan to 
continue the search for a good reflectance function. Someday 
we may even do it. 

3.5 Computation of Sun Position C^) Prom Ephemeris Data 

In order to compute the sun's position angle, a , we must know: 

1. where the sun is, 

2. where the moon is, 

3. where the sensor is, 

all with respect to the center of the earth. To accomplish 
this, we introduce another coordinate system, E , which is 

tied to the celestial sphere and has its origin at the center 

c c 

of the earth. E^ points to the north celestial pole and E£ 

is aligned to the First Point of Aries. Thus, the celestial 

sphere right ascension variable has 0^ aligned to E^ and 6^ 
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Figure 3-8 

Relationship o£ Coordinate System 
E c to the Celestial Sphere 




aligned to E^; ref. Figure 3-8. The earth rotates about 

C C 

and Greenwich lies in the positive - E£ plane at 0 U.T. 

Suppose a point on the celestial sphere has a right ascension 
and declination of p and 6 respectively. The unit vector 
aligned to this point will have the form 


s6 

cS cp 
c<5sp 


3.34 


A point on the earth having an east longitude of X and a lati- 
tude of <j> can be represented by 


" s<|> - 
c<j)cX 
C(j>sX 


3.35 


in an earth-based coordinate system. To compute its coordinates 
in E we must account for the -fact that the vector is rotating 
about E-^ at an angular rate oj and has been doing so for some 
time t. The transformation is performed as follows 
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”l 

0 

0 


_ _ 
scj> 

0 

cut 

-sut 


ccf>cA. 

0 

swt 

cut 

c<f>sAj 


3.36 


Note that this 3x3 matrix rotates vectors about the 1-axis 
in a fixed coordinate frame . 

Now we can proceed. Look at Figure 3-9. It depicts the per- 
tinent elements of the analysis; earth, moon, lunar sensor, 
and sun. The view is from above the plane defined by the 
center of the moon, the center of the sun, and the lunar sensor. 
Note that this plane is none other than the 2-3 plane used in 
earlier sections; e.g., look at Figure 3-3. The vectors G, 

L, and H extend from the center of the earth to the sensor, 
the lunar center, and sun’s center respectively. As such, 
they do not necessarily lie in the plane defined by D and S* 

(our familiar 2-3 plane) . One point of clarification; the vec- 
tor S*, appearing in Figure 3-9, differs from the vector S, 
Figure 3-3, only in magnitude. S is a unit vector. 

Our problem is to find o and we can do this if we know D and 
S* for then 
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0 


-1 


D-S* 


3. 37 


= cos 


D | Is* 


To calculate D, we first need G and L. If h is the distance 
from the center of the earth to the sensor, then 



”l 

0 

0 

s<f> 

G = h 

0 

cut 

“SUt 

C(j>cX 


0 

sut 

cut 

1 C(j)sX 


If d is the earth-moon distance and p and 6 are the moon's 
m m 

right ascension and declination, then 


L 



r 


s<$ 


c<$ 


m 


cS 


m 


m 

cp 

sp 


J 

UlJ 


From Figure 3-9 we see that 


D = G - L 


3. 38 


3.39 


3.40 
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Next, to compute S*, we need H in addition to L. Let d g be 
the earth-sun distance; p g and are the sun's right ascension 
and declination. Then 


H 


d 


s 




3.41 


Again, from Figure 3-9, we see that 


S* = H - L 


3.42 


Thus, given the sensor's geographic position, the sun and moon's 
celestial coordinates, and the time, equations 3.38 through 
3.42 can be used to compute D and S* which in turn, via equa- 
tion 3.37, produce a. The sun and moon data can be obtained 
from an ephemeris and the particulars of doing so are discussed 
in Section 4. 

3.6 Computation of Sensor Offset Angle 

Once we have computed all the pertinent quantities developed 
in the preceding sections we will know where the sensor will 
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point... in terms of a point on the lunar equator which is dis- 
placed an angle 0* from the 2-axis. All that is left is to 
convert this into an angular offset of the sensor about its 
yaw axis . 

Consider Figure 3-10. It is a top view of the 2-3 plane and 

shows the moon, the vector D (ref. Section 3.S), and the edge- 

view of the sensor's yaw-null plane. The sensor is nulled and, 

as such, its yaw-null plane intersects the lunar equator at a 

* 

point which is displaced an angle 0 O from the 2 -axis. 



Figure 3-10 

Geometry for Calculating Sensor Offset Angles 
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Recall; D is the vector from the center of the moon to the 
sensor and has magnitude d , r is the radius of the moon. 

From Figure 3-10 we have 

r s0o = z = (d - r c0 o ) tan e 3.43 

m 


or 


e 


tan 



r s9o 

£ 

- r c9 o 


3.44 


■■This quantity, e, is the angle the sensor will be displaced 
from the center of the lunar disk when the sensor is nulled. 
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Section IV 


COMPUTER PROGRAM DESCRIPTIONS 
4 . 0 INTRODUCTION 

The purpose of the computer programs written for the lunar 
radiometric center problem was, in a nutshell, to evaluate 
the predicted vehicle sensor offset angle when pointed pre- 
cisely at the moon's center of brightness rather than its 
physical center. The solution to the problem involved, pri- 
marily, determination of required variables to evaluate the 
integral of the power density function for a specified time 
and location of launch, and from the results of that integra- 
tion, to calculate the radiometric center of the moon and 
thereby determine the predicted offset angle. 

These calculations were accomplished in one mainline program, 
assisted by three subroutines, written in FORTRAN IV language 
for the BBRC engineering computer, a General Automation 18/30. 
Descriptions of the general logic flow of each of the programs 
constitute Section 4.2 and are preceeded by some explanation 
of features common to all of the programs in Section 4.1. 
Section 4.3 provides listings and examples of output for each 
of the programs. 
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4.1 


General Notes 



4.1.1 Accuracy of Calculations and 18/30 Computer Precision 

Because of the relatively complex nature of and the many itera- 
tive calculations required by the lunar radiometric center pro- 
blem, the computer programs were very susceptible to error 
propagation. The use of the extended precision feature signi- 
ficantly reduced error build-up problems, but a substantial 
effort was also made through coding and numerical analysis 
techniques to avoid inaccurate results due to error pitfalls. 
This effort was particularly applied in the task of incre- 
menting loop variables when the increment values were frac- 
tions which could not be exactly represented in binary. In 
this case simple accumulation of the variable value (x = x + 

Ax) also results in accumulated error, but the problem was 
avoided completely by always returning the variable to a base 
value before adjusting it for the next iteration. 

All the programs performing the lunar radiometric center cal- 
culations were run in G.A. 18/30 Extended Precision, which 
utilizes 3 16-bit words to represent real variables. Mention 
is made of this fact to differentiate it from the more commonly 
known "double precision" which uses four-word floating point 
representation. The 18/30 is not provided with the double 
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precision feature, but the use of extended precision yields 
between nine and ten significant digits of accuracy which was 
adequate for these calculations. 

4.1.2 Computer Program Execution Time 

Again because of the iterative nature of these calculations, 
it was necessary to reduce, as much as possible, the computer 
run time required to complete the computations. The problem 
was aggrevated by the fact that the 18/30 is not supplied 
with floating point hardware, but rather must rely entirely 
on software subroutines to perform all floating point opera- 
tions. The use of extended precision intensifies the run 
time problem even further. The resulting calls to floating 
point subroutines, especially in the long integration loops, 

greatly increased run time requirements. As a result, exten- 

\ 

sive efforts were made in the programming of the lunar radio- 
metric center calculations to simplify equations in loops in 
such a way as to reduce the number of redundant subroutine 
calls and thereby decrease program execution time. Another 
solution was to integrate only over half of symmetrical regions. 
Because the majority of program calculations were involved with 
the integration of the power density function, use of this 
technique reduced execution time by a factor of two. (See 
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Subroutine BIGNT description for further discussion of this 
integration technique) . 

4.2 Program Descriptions 

4.2.1 Mainline Program, JAKE 

The primary function of the mainline program, named JAKE, was 
to initialize parameters required by the power density func- 
tion integration subroutines. Upon return from those sub- 
routines, JAKE operated on the integration results to calcu- 
late the anticipated sensor offset angle. 

4. 2. 1.1 Parameter Initialization 

Five parameters were required as input to the integration sub- 
routines: the lunar phase angle, earth-moon distance for the 

specified launch time, radius of the moon, and increment values 
for the two variables of integration. The lunar radius was con- 
stant and so was initialized as data. The two integration 
variable increments values were constant throughout the cal- 
culations and were read from card input. The other two para- 
meters were calculated using ephemeris data for the launch time 
and latitude and longitude of the launch site. 

a. Calculation of Earth-Moon Distance 

The American Ephemeris provides the necessary equa- 
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tion for calculation of the true geocentric distance 
for given time of day: 


D 


(a Q + a-jP + a 2 P 2 + a 3 P 3 ) 6378 . 16 


4.1 


where 


p is a number between zero and one representing the 
fractional part of the half-day in which the GMT lies. 

a , a^ & 2 > a 3 are polynomial coefficients for true 
geocentric distance, tabulated in the ephemeris for 
each half day of the year (expressed in units of 
Earth's equatorial radius). 

The constant 6378.16 is the earth's equatorial radius 
in kilometers and is used to convert the earth-moon 
distance into kilometers for consistency with the 
other calculations . 

b. Calculation of the Lunar Phase Angle 

The vector equations required to calculate the 
lunar phase angle were described in Section 3.5 of 
this report. Parameters required for these calcu- 
lations are: 
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Sidereal time for 0 hours universal time of the launch 
date; right ascension, declination, and geocentric dis- 
tance of the sun; right ascension, declination, and geo- 
centric distance of the moon; and latitude and longitude 
on the launch site. 

Latitude and longitude of the launch site and sidereal 
time were read as card input requiring no further adjust- 
ment and the calculation of lunar geocentric distance was 
described above. Right ascension, declination, and geo- 
centric distance of the sun were obtained from the Ameri- 
can Ephemeris for 0 hours Ephemeris Time of the days pro- 
ceeding and succeeding the launch, and right ascension 
and declination of the moon were obtained from the same 
source for the hours of Ephemeris Time preceeding and suc- 
ceeding the launch. These values then needed to be ad- 
justed to reflect the exact time of the launch, which was 
accomplished via the following equation for general linear 
interpolation: 

y = ( x Cy 2 - + x 2 y i ~ x iy2 )/(x 2 * x i> 

where y = the interpolated right ascension, declination 
or distance 
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x = the Universal Time of the launch (hours and 

minutes for solar variables, minutes for lunar 
variables) 

and x 2 = 0 and 24 hours, respectively for solar variables, 
and 0 and 60 minutes, respectively, for lunar 
variables . 

y^ and y 2 = right ascensions, declinations, or distances 
as obtained from the ephemeris for the times 
bracketing the time of launch. 

See Figure 4-1 for a graphic illustration of lunar inter- 
polation. 

Once these interpolated values were computed, calculation 
of the lunar phase angle was accomplished as described in 
Section 3.5. 

4.2.1. 2 Calculation of the Sensor Offset Angle 

The calculation of the sensor offset angle is described in Sec- 
tion 3.6 of this report, and of necessity, some of the variable 
names mentioned there will also be used here. In order to 
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CO 


y 



TO INTERPOLATE , ASSUMING Y = f { x ) IS LINEAR: 
Y=(X(Y 2 -Y|}+X 2 Y| — X| Y 2 )/ (x 2 -x 1 ) 



Figure 4-1 

Generalized Linear Interpolation 


determine the Sensor Error Angle, it was first necessary to 

* 

determine the value for 0o which would correspond to the lunar 
radiometric center, as shown by the following ratio: 

B C9 o ) _ i t -> 

BO/2) ” 1/2 

The power density function integration subroutine returns to 

the mainline an array of tabulated integration values (B(6 0 ) 's) 

corresponding to the range of G 0 angles. Dividing the end 

point of this array (B(tt/ 2)) by two yields the center of bright- 
* 

ness (B(0 O ) value. By finding the two elements in the array 
which bracketed that point, it was then only necessary to inter- 
polate between the two corresponding 0o values to determine 
* 

the required 0 O angle. This interpolation was accomplished 
using the linear equation described above. Having found 0 O , 
the Sensor Offset Angle was calculated as shown in Section 3.6. 

4.2.2 Subroutine BIGNT 

The purpose of subroutine BIGNT x^as to perform the required 
integration of the power density function, given the necessary 
parameters provided by mainline program JAKE. BIGNT 1 s first 
task was to determine the limits of the integration variables, 
a and 0 O . The criteria for this determination are described in 
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Section 3.3 of this report. Initially, the incremental values 
of integration, Aa and A6 0 , were fixed at 1°, and since the 
upper and lower bounds of 6o were always integral amounts of 
degrees, the 1° restriction on A6o was sufficient. The inte- 
gration limits of a, however, often included fractional parts 
of degrees, so the Aa value was subsequently adjusted slightly 
to fit exactly within the range of integration. 

It is also necessary to note that the initial values for the 
variables of integration were not set exactly to their respective 
lower limits, but rather were set a half increment above those 
values (i.e., lower limit + Aa/2 or +A0 o /2). The reason for 
this initialization was twofold: 1) Since only half of the 

full symmetrical region was being integrated. CSee Section 
3.3), the variables of integration were set up a half step to 
fully cover the illuminated region of the moon and to avoid 
duplication of the region's mid-area; and 2) the most meaning- 
ful result of the point-wise integration could be obtained by 
evaluating the integrand at the mid-point, rather than an end- 
point, of the differential area. See Figure 4-2 for an, illu- 
stration of this method of area mid-point integrand evaluation. 

Having accomplished all of this initialization, BIGNT was then 
free to attack the problem at hand, the evaluation and integra- 
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Figure 4-2 

Illustration of Mid-Point Integrand Evaluation 




tion of the lunar power density function. The task was accom- 
plished in two nested DO loops, the outer loop being concerned 
with 0o, and the inner loop with a. For each value of 0 O , the 
power density function was evaluated for every value in the 
range of a and the resultant function was then integrated over 
a. The actual integration was performed by an 18/30 Scientific 
Subroutine Library program, QSF, which uses a combination of 
Simpson's Rule and Newton's 3/8 Rule. QSF operates on a tabu- 
lated array of function values, and returns a tabulated array 
of step-wise integration values, the endpoint of which is the 
final result of the integration. 

M 

As each a loop is completed for each 0o value, the endpoint 
of the resulting integration vector was used to form an ele- 
ment of the array which, when the 0 O loop was complete, would 
again be integrated by the QSF subroutine, this time over 0o. 
The final array of integration values, thus obtained by BIGNT, 
was then returned to the mainline program along with an integer 
variable indicating the number of elements in the array. Hav- 
ing completed its appointed rounds , BIGNT then goes to sleep 
awaiting another call from the mainline program and a new set 
of parameters. 
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4.2.3 


Function Subprograms FLUX and RDFCN 


The actual evaluation of the power density function at a single ' 
point was performed by two small function subprograms , FLUX 
and RDFCN. Subprogram FLUX was responsible for evaluating the 
entire power density function, and it calls RDFCN for defini- 
tion of the lunar reflectance function. During the development 
of the math model, when experimentation with Lambertian, Lommel- 
Seeliger, and other radiance functions was going on, it was 
most convenient to have the function evaluation and reflectance 
function itself separated from the larger programs so that small 
changes to the equations could be easily and quickly incorporated. 
Since that experimentation continued up to the very last calcu- 
lations, RDFCN and FLUX remained in existance and were never 
merged with BIGNT. 

4.2.4 Additional Programs 

Using BIGNT, FLUX and RDFCN as a calculating base for the lunar 
radiometric center problem, several other mainline programs , 
besides JAKE, were written. Since these programs did not con- 
stitute part of the final product, but only contributed to its 
formulation, they will not be described in any detail here. 

They should, however, be mentioned because of the functions 
they provided as constructive tools. 




Primarily, the other mainline programs differed from JAKE in 
that they did not actually calculate the lunar phase angle and 
earth-moon distance parameters required for the integration of 
the power density function. Instead, they operated on constant 
values for these variables. Neither were those programs con- 
cerned with the determination of the sensor offset angle, but 
rather centered around the calculation of power density ratios 
for different lunar phase conditions. No consideration was given 
to an actual launch situation, but instead hypothetical, and 
therefore somewhat predictable, circumstances were assumed. 

In this sense, these programs provided a valuable tool for 
evaluating the performance of different lunar reflectance func- 
tions when comparing, for instance, full to half-moon bright- 
nesses. These programs also provided a wealth of plotted out- 
put which served as a visual aid in the trouble-shooting of 
both the programs and the math model. So although these pro- 
grams did not constitute a part of the final calculations re- 
quired by the lunar radiometric center problem, their contri- 
bution to the end item cannot be overlooked. The listing for 
mainline program HERB is included in the next section to serve 
as an example of the type of function provided by these inter- 
mediate mainline programs. 
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Program Listings and Sample Output 

Included in this section axe listings and sample output, where 
applicable, for each of the programs and subroutines which con- 
tributed to the solution of the lunar radiometric center pro- 
blem. The order of appearance is as follows: 

1. Mainline Program JAKE 

Sample teletype output 
Sample printer output 

2. Subroutine BIGNT 

Sample teletype output 

3. Function Subprograms FLUX and RDFCN 

4. Intermediate Mainline Program HERB 

Sample plotted output. 


4-15 




LISTING OF 

MAINLINE PROGRAM JAKE 
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PROGRAM JAKE* EXERCISES LUNAR BRIGHTNESS FUNCTIONS FOR ONE SIGMA 


C-ERRS. . .STNO.C. 


FORTRAN SOURCE statements 


identfcn 


♦♦COMPILER MESSAGES^ 


CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 
C PROGRAM NAME - JAKE- ~ C 


DETERMINES LGNaR PHASE ANGLE, SIGMA, FROM LUNaR 
and solar RIGHT AS CEN SIONS, DECLINA TIONS, Ea RTH 

distances, etc. "for th’e specified" time and loca- 
tion of THE LaUNs.H SIGMft IS THEN O UTPU T TO SUB 
ROUTINE 'BIGNT't WHERE THE LUNaR BRIGHTNESS FUNC 
TION IS EVALUATED. UPON RETURN TO »JArE», THE 


LUNaR CENTER OF BRIGHTNESS angle, THETa-0 IS de- c 

terminedi from Which the Jr acker o ffset a ngle* e P- c 

" SILON is" CALCULATED, C 

_ c 

SUBROUTINES REQUIRED m *816 NT » ,* BRIGHTNESS FUNCTION TnTeg(<aTOR c 

• rot * , three-dimensional coordinate sys- 


tem ROTATION SUBROUTINE 
•DMS2R* AND *HMS2R*« CONVERSION ROUTINES 
FOR EPHEMERIS HOURS AND DEGREES TO 
. _ _RAD_IANS 

CARO INPUT - LAUNCH DESCRIPTION HEADER CA«D 


RIGHT ASCENSIONS, DECLINATIONS AND EaRTH OIS 
tances OF SUN FOR DOY AND DOY+1 

SIDEREAL TIME_F0R 0 HOURS ON DOY 

RIGHT' ASCENSIONS, “declinations and EaRTH DIS 
tances of moon for hours and hours+i 


MOON-EARTH DISTANCE POLYNOMIAL coefficients 
latitude and longitude of LAUNCH site 




BALCbROTHERS research CORPORATI 
BOULDER, COLORADO SEPTEMBER « 1974 


cccccccccclccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 


c 

COMMON SIGMA, RAD, DM 


REAL INTRP * LAMDA 

INTEGER DOY 


C 

DATA AU/1.496E6/, SF/l.OE-6/ ’ 


DATA OMEGA/15 *0410 7/ 

DATA PI/3. l*nS926Sif/ 


DIMENSION ISTUF<36> *aC+> *BM(3>,bS<3> »8V<3> iBD(3) ,E0<3>*BE<3) 


DIMENSION 6(180) 


R2D(X) = X*180.0/PI 
D2R(X) ts X*PI/180.0 




ROUND ( X ) 


IFIX(X*XO0,O + 0.5 ) /100 . 0 
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program Jake, exercises lunar brightness functions for one sigma 


,,,STNO.C..... FORTRAN SOUR C__E S TAT E M E N T S 

CCCCCLCLCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCIXCCCLLCCCC 

c c 

C DO SOME INITIALIZING “ 


ioentfcn 


**C0MHILE« MESSAGES** 


call TIMEX<0,LL,KK) 
RAO = 17358.0 


READ(2,4000,ENU=500) ISTUF 
REAO( 2 , 4001 ) OOY » IH, M 
WRITE< 1,8000 ) 

WRITE! 5,8001) 

WRITE ( 5 , 4007 ) ISTUF 
WRITt(l,4000) ISTUF 
WRITE(5,<tlOO) OOY* IH, M 


UT = IH + M/60.0 
P = UT/12.0 
IF(IH ,GT» 12) P 


P-1.0 


CCCCCCCCCLLCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCLCCCCLCCCCC 

C DO SUN’S R«a*» 0£C»* A n D DISTaNcE INTERPOLATIONS, . .READ THETO C 


REAU<2,40C2) IUAY,IRHl,IRMl,RSl,IOHl,IUMl,USl,SOSl 

REAU <2, 4002) L0AY«IRH2,IRM2,RS2,IDH2,IUM2,US2,S0S2 

IF(IUAY .EQ, DOY .and. LDAY ,EQ. OOY+l) go to 2 

WRITE(S,4005) 

PAUSE 1234 

GO TO. 1 - -- 

READ ( 2,4009 ) IT_H,ITM,TS_ „ _ — --- 

THEfO = (ITH + (ITM + tS/60,0)/60,0)*15,0 

WR IT ET57420T) ~ " ~ 

WRITE ( 5, 4300 ) lRHl,lRMl,RSl,IDHl_,.IpMl,DSl,SDSl _ , 

WRITE ( 5 , 4300 ) IRH2,IRM2*RS2,IDH2,IUM2,DS2,SDS2 

WRITE 1 5 , 4400 ) ITH»ITM,TS _ , ... - - 

SR A = INTRP(Q t 0,24.0, HMS2R( I Rhl,IRMl, RSI), HMS2R(IRH2,IRM2,RS2),UTj 
SDEC = INTRP{0.0,24.,DMS2R<lDHl*IDMl,DSl|»DMS2R( IDH2 , IDM2,DS2) ,UT) 
OS = I NTRPtO, 0,24.0, SDS1,SDS2,UT) . 

DS = DS* AU c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

DO INTERPOLATIONS FOR MOON'S R.a. AND DEC., AND COMPUTE DISTANCE C 
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PROGRAM JAKE, EXERCISES LUNAR BRIGHTNESS FUNCTIONS For oNE SIGMA 


C- ERft S. 4 .STN 0 .C... v « F 0 r_T K A N S 0 y r C E ST A tEM E NTS 


IDENtFCN 


**COHpICEf< MESSAGES** 


C FRON MOON TO EARTH USING A-COEFFICIENTS A w O P C 

C C 

ccccccccccuccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

3 . READ ( 2.4QQ2 ) IN, IRHl , IRMl, RSj., IDHI^IQMj., D$1 

R€AD<2tit00S> LF! t JRH2* IRM2* RS2* IDH2* IDM2* 0 S 2 
C 

IF(IM .EQ. IH .AND. LM ,EQ. IH+1) GO TO 4 
WRITE (5.4005) 

PAUSE 2345 

GO TO 3 _ . . . 

READt 2.4004 ) A 
WRITE ( 5.4500) 

WRITE 15*4600 ) IRHl* IRMl * RSl * iDHl , IOMl * DSl 

WRITE (5, 4600) IrH2,Irm2«KS2*IDH2,IDm2.DS2 

WRITE ( 5*4700) A 

RAM = INTRP<0.0,60.0,hmS2R(IRh1«IRM1*RSI| *HP1S2R<IRH2,IRM2*RS2>,FL0 
1 AT ( M ) ) 

D_EcM = INTRP ( 0 . 0 ,60,0 .DMS2R ( luHl t IDMl 1 uSl > , £)MS2 r ( I0H2 , Iul"|2 * US2 ) ,FU 
10 AT t M) ) ' " 


C 

4 

C 


10 

c 


DO 10 0 = 2,4 
A< 0 ) = A < 0 )*SF 
CONTINUE 


DM = UfAC**‘)*P + A(3'))*p + A < 2 ) ) *P + A(l) )*63TF7l6 
C 

WRITE ( 5 1 4006 ) SRft, SDEC, DS, RAM, DECM, DM 

C C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

'§ cTrcULATE'-rEcrTrRs'M-, -ST “£0 ,~E 7-V-JIW0 §- 

C __ c 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


BM ( 1 ) = SIN(DECM)*DM 

BM< 2 > = COS(OECM)*COS(RAM)*DM 

bmT3'1~=~ 'Cos ( de"Cm ) *sin (ramTSDm 


BS ( 1 ) = SIN ( SDEC ) * 0 S 

BS ( 2 ) = COS<SOEC)*COS<SRA)*DS 

BS ( 3 ) = COSlSDtC)*SIN(SRA)*DS 


READ ( 2 , 4000 ) LAMDA, PHI 
EO ( 1 ) = SIM< D 2 R t PHI ) ) 

EO ( 2 ) = C 0 S(D 2 K(pHI) )*COS( 02 R(LAMOA) ) 
EO < 3 ) = COS(D 2 R(PHI) >*SiN{D 2 R(LAnbA) ) 


c 


THETA = D 2 R( 0 MtGA*UT + THETO) 
CALL ROTI-THETA, 1 , EO* BE, IER) 
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PROGRAM JaK£« exercises lunar brightness functions for one sigma 


C-ERRS...STNO.C FORTRAN 


SOURCE STATEMENTS^, . , ... IOE NTFCN ♦♦ COMPILER fltbSAb ES** 


THET = OTHET/2.0 

IFISIGMA .LE. 90.0) THET = THET + SlGflA - 90.0 
THETA = THET 


-F* 

i 

to 


c 


11 

GO 

C 

22 

C 

23 

C 

70 


C 


C 

c 

500 


C 


DO 60 1=1, J 

IF ( 6 ( I ) - HALF} 11,22,23 
A1 = THETA 

THETA = THET + I*DTHET 
CONTINUE 

AO = THETA 
GO TO 70 


A2 = THETA 

AO = (HALF*< A2-A1) + B ( I ) *A1 


B < 1-1 ) #A2 ) / ( B ( I ) 


B(I-l) ) 


EPS = ATAN<RAD*SINf02 R< a< 5 ) )/(OM - RAD*COS < 02R ( ftO ) ) )) 
EPS = R2D ( EPS ) *60 , 0 
EPSI = ROUND (EPS ) 


WRITE ( 5 , 6200 ) AO 
WRITE ( 5 , 6300 ) EPS 
WRITE ( 1 , 6900 ) DOY, IH, M 
WR ITE( 1,6600) PHI, LAMOA 
WRITE? 1 , 3000 ) SIGMA 
WRITE ( 1 , 600 0 ) EPSl 


GO TO 5 . . 

CALL TIMEX(1,LL,KK> 

XK = KK/100,0 

WRITE ( 5, 6100 ) LL. XK 


_£ALL EXIT _ . „ 

"C” 

C - -? 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCO 

C ... _ . . 

C * FORMAT STATEMENTS c 

c £ 

CCCCCLCCCCLCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


c 


c 


C 

c ...... 

3000 FORMaT(//’LUNAR phase angle = * ,F8.4»* DEGREES*) 

C 

'4000 FORMAT (36A2 ) 

C 

4001 F0RMAT(I3,2(1X,I2> ) 

C 

4002 FORMAT < 13, IX, 2 (14, 13, IX, F6. 3) ,F10,7) 

C — 

"4003 FORMAT ( *OSCREW Up ON DOY FOR SUN DATA. ..TRY AGAIN.’) 

C 

4004 FORMAT ( 4F10 • 0 ) 



zz- 
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PROGRAM JAKE, EXERCISES LUnAr BRIGHTNESS FUNCTIONS FoR ONE SIGMA 


c-Errs...stno* c »**** fortran source state mi e n_t _s 


IDEnTFCn 


♦♦Compiler MEs>sAGts*^ 


c 

4005 FORMAT( *oSCREW UP ON HOURS FOR MOON OATA...TRY AGAIN*) 

o 

4006 FORMAT {//' 0*,*SUN'/* '*5X,*R.A* = '*£14.7,* DEC. = '»El4,7»* AND E 
IArtH DISTANCE = «,F 14 . 7 /* 0 *j*MOON*/* « ,Bx, 'R. A. =_*,E14.7_,' DEC. = 

' 2 »',E14,7,' AND EARTH DISTANCE = ',Eli).7/) 

C 

4007 FORMAT { » 1 ’ , 36A2) 

C 

4003 FORMAT12F10.5) 

C 

4009 F0RMAT<I4,I3,1X,F6.3) 

Q 

4100 FORMAT! ’ 0 ** 'qOY ’ i3,5X*' TIME s ' ♦ 12 * * H * * 2X , l2 , * m * UNIVERSAL TIME 
2 ’ ) 

C 

4200 FORMAT(//*o* , ’EPHEMERIS DATA, SUN*) 

C 

4300 FORMAT! *0* »5X* 'RIGHT ASCENSION s * , 14 » 13 , IX * F^ . 3 , 2X * ' * DECLINATION 
1 = * , I4,I3,lX«F6.3, • AND DISTANCE = *,F10,5) 

C 

4400 FORMAT !//' 0 * , * EPHEMERIS THETA-0 ANGLE = • * 14 » I3,1X* F6»3) 

4500 FORMAT !//*0 * , 'EPHEMERIS OATa, MOON*')*" 

c 

4600 FORMaT(»0'*Sx*'RIGHT ASCENSION = »*i4,i3,1X,F6,3,2X,' AND DECLINaT 
1ION = »,I4,I3,1X.F6.3) 

Q 

4700 FORMAT! • 0 * ♦ 5X_* * M OON Q ISTaNCE P OLYNOMIAL C OEFFICIENTS = »,F 9.6* 3!2X 
1,F10.1)) 

0 

5000 FORMaT!//' 0* ,5X, *M VECTOR s * * El4 . 7 , 5X* • S VECTOR = » ,e14,7* SX* *EO 
XVECTOR = ♦ ,E14.7) . - 

C 

_ 50 0_1_ F0RM AT( 2( *0»,16X. E14 .7, 16X*E1‘».7*17X *E14.7/ )) 

c 

5002 F0RMAT(//*0* ,5X, •£ VECTOR = *, E14 . 7 , 5X ,* V RECTOR = *_, £14 %1 * 6X,_* D V 
1ECTOR = *,E14.7) 

5003 ‘FORMAT!//*!)*, 'SIGMA = *,E14",7) ~~ 

5TJTJ4 — FORMa T~r/ / r tHT ' L ' fl TTTDDTr~TfTD~CPN'GTTTUPE'~Oy LfiUNUF'SITg = * ,£14", TV~ 

1 AND * , E14.7 , * * RESPECTIVELY' > 

C 

6000 FORMAT! 'TRACKER OFFSET ANGLE = **F5. 2, • .ARC-MINUTES’ > . 

C 

6100 FORMAT(//'0* , 'EXECUTION TIME = « , 13 . IX , F 5 . 2 ) 

C " 

6200 FORMAT! 'O' * 'THETO* CENTER OF INTENSITY ANGLE = *,El_4.7) 

C 

6300 FORMAT! *0' * 'EPSILON* LOS ANGLE = **E14.7,* ARC-MINUTES') 

C 

6600 FORMAT!//'LATITUDE OF THE LAUNCH SITE = »*F7,3,' DEGREES' /'LONGIT 

luce U K TH E T W ITC H" SUE = ' « F a , S , * D E GREE ' S ' ' ) ’ 

c 

6900 FORMAT!//’DOY = '*13,' AND LAUNCH TIME = ' 1 12* ’ H* * IX, 12 * *Mt UNIVER 
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program Jake* exercises lunar brightness functions for one sigma 


C-ERRS...STNO.C FORTRAN SOURC E ST A T t M E N_ T S 


IOENTFCN 


♦♦COMPILE* MES&AGtS** 


1SAL TIME' ) 

C 

6000 F0RMAT</////> 

C 

8001 FORMAT(’I') 

c 

c 

END 

VARIABLE ALLOCATIONS 

blank common block 

SIGMa(R*6 c>=7FFO RaD(R*6 C>=7FFA 


EQUIVALENCES 
LAMDA ( R*6 
IH { 1*2 
IRHl ( 1*2 

OSl (p*6 

RS2 ( R*6 
ITH ( 1*2 
sDeC<r*6 
RAM(R*6 
THETA (R*6 
PR00(R*6 

HACF'j R*6 
A2<R*6 
A<f?*6 
E0(R*6 


a INTERNAL 

)=0g00 

>=0209 

>=0212 

>=0219 

>=0222 

>=0220 

>=0236 

>=0243 

>=0250 

)=025D 

>-0260 

>=0279 

>=02DA»02D1 

>=0307-0301 


STATEMENT ALLOCATIONS 
3000=”05S2 "4000=U576 

4O0 9 =O6oC 4100=0611 

5002=06FB 5oo3=u7lB 
6O0l=07FE R2D=0801 

10=0A5A 30=0B36 


variables 

DOY ( 1*2 
M(I*2 
IRM1 ( I*2 
S0S1(K*6 
IDH2 < 1*2 
UM( 1*2 
DS(R*6 
DECM ( R*6 
OMEGA ( R*6 
OeNOM { R*6 
THET ( R*6 
EPS { R*6 

»M(R*6 

BE ( R*6 


400 1 = 0^5 79 _ 
42QO=o62E 
5004=0726 
D2R=0B0E 
4o=qB72 


O 2 O 3 

020A 

0213 

021C 

0225 

022E 

023B 

0246 

0253 

0260 

026F 

= 0?7C 

=02^3-0200 
=0310-030 A 


DM<R*6 C>=7fF7 


PI ( R*6 
UT(R*6 
RSl ( R*6 
LDAY( 1*2 

rDM2(I*2 
TS ( R*6 
Au { R*6 
U( 1*2 
IER< 1*2 
0AlPH(R*6 
T (T*2 
EPSI < R*6 
Bs(R*6 
B ( R*6 


4*0 0 2=0'57f JTOO 3^o5W ' 

43oO=0 6 3E 4400=0668 

6o0O=074F 61 00=0 765 

INTRPsOBlB ROUND=064 c 


11=qC21 


6o=oC2E 


sOgOH 
= 0208 
= 0214 

= 021F 
= 0226 
= 022F 
= 023E 
= 0249 

=0256 
= 0263 
= 0272 
= 027F 

=02EC-02E6 

=0520-0313 


"WtFfrCfDA 0 
45oo=067F 
6200=6777 
5=0883 
22=0 c 37 


LL 1 1*2 

>=0207 

KM 1*2 

>=0208 

P ( R*6 

)=020£ 

IUA y { 1*2 

>=0211 

IDUl ( 1*2 

>=0217 

I0Ml t I*2 

>=02lS 

IRFI2 ( 1*2 

>=0220 

IKM2(1*2 

>=0221 

□S2(R*6 

>=0227 

SDS2{R*H 

|=022A 

THETO < R*6 

>=0232 

bRA(K*fa 

>=0235 

IM( 1*2 

>=0241 

LM ( 1*2 

>=0242 

SF ( R*6 

>=U24A 

PHI(K*6 

) =0240 

DENitR*6 

>=0257 

DLN2(R*6 

>=025A 

DTHET ( R*6 

>=0266 

BBl ( k*6 

)=026g 

Al(R*6 

T=W3 ’ 

Su7R*5 " 

> soars 

XK ( R*6 

>=0262 

ISTUF t 1*2 

)=O2U0-02AD 

Bv<R*6 

)=02F5-02EF 

8D{k*6 

)=02FE-02F8 


’" i t00’5'=0^ A 3~ 
4600=068F 
63qo=o70£ 
1=08 d 6 
2 A = 0 C 3D 


4700=g6Bl 

6800=67115 

2=0916 

7o=6‘CbF 


5000=Q6CF 
6900=07Qb 
3=09A6 
500 = 0 t - Ctt 


“'4 0trB'="o?O 9 - 

5l)0l = 0&EF 
8uo6=07T8 
4=09E2 


FEATURES SUPPORTED 
ONE WORo INTEGERS 
“ EXTENDED PRECISION 
ORIGIN , _ . 

IOCS- 

PLOTTE R 

1403 PRINTER 

TYPEWRITER 

CARD 

CALLED SUBPROGRAMS 

EIFIX TIM E X HMS2R pMS2R EFLOT _ESI_N 


EsUSX Empy 'EfipYx 

SWRT ScOMP SFIO 

LRGT LRLE LRLQ 

REAL CONSTANTS 

.180000000E 03=0530 
7120000000L 0 2=0 53F" 
.63781600QE q4=054E 


SIOAI SIOAF SiOF 


• i ooouoo^E 

IhoOOOOOUOE 01=0551 


ROT esqrt eacos bIGNT LaTa n INITS EADLJ esub 

"TsYoS — CsBr edvh nre I GFrTTZ — cmz — vchri sW 

sugsc prnz pause snr subin seof float xsub 


.500000000E 00=0536 

"^L^OFOOOUOE "02=0545 

,9UOUOOOOOE 02=0554 


.173800000E 04=0539 


.60OUOUOUOE U2 s 053c 
r2TP013WruTfe~02=0'59B- 



LAUNCH AT KAUAI, HAWAII, 3 N0V. 1974, I OH 38M U.T. 

D0Y = 307 AND LAUNCH TIME = 1 OH 38M, UNIVERSAL TIME 

LATITUDE 0F THE LAUNCH SITE = 22*065 DEGREES 

LONGITUDE 0F THE LAUNCH SITE = -159*781 DEGREES 

LUNAR PHASE ANGLE = 42*1534 DEGREES 

TRACKER OFFSET ANGLE = 4*88 ARC-MINUTES 


EXAMPLE OF FINAL TELETYPE OUTPUT 
FROM MAINLINE PROGRAM JAKE 


4-24 




EXAMPLE OF PRINTED OUTPUT 
FROM MAINLINE PROGRAM JAKE 


4-25 
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LAUNCH AT KAUAI, HAWAII, 3 NOV. 1974. 10H 381*1 U.T. 
DOY307 TIME = I0H 38M» UNIVERSAL TIME 


EPHEMERIS DATA* SUN 

RIGHT ASCENSION - 14- 31 4.990 * UECLINaTION = -14 52 24.700 AND DISTANCE = 0.99203 

RTghT ASCENSION'S *14 35 'i ~37 cT ", DETCTRaTT AIT'S rj.1T 1.1 lo.'800 AND "DfSTANC£”= oTVS 178 


EPHEMERIS THETA-o ANGLE = 2 47 29.308 


EPHEMERIS DATA. MOON 

RIGHT ASCENSION = 5 27 36.739 AND DECLINATION = 22 8 15.320 

RIGHT ASCENSION - 5 30 5.930 AND OtCLlTTATION = 22 7 2,450 ~ 

MOON DISTANCE POLYNOMIAL COEFFICIENTS r 58.992645 -196552.0 10632.0 -39.0 


SUN 

R.A. = 

0,3Bd8433E 

01 

DEC. = - 0 . 2620108 E 00 

AND 

EARTH 

DISTANCE a 

0.1483921E 

09 

MOON 

R* A* " 

0.1436349c 

01 

DEC. = 0.38615001 00 

AND 

earth 

DISTANCE = 

0.3752067E 

06 



LATITUDE ANO LONGITTUoe OF LaUNcH SITE = 0.2206553E 02_ aNQ -0.1597818E 03, RESPECTIVELY 


M VECTOR = 0.1413121E 06 

0.465901lE_p5_ 
0 . 3444420E Qfe 


S VECTOR = - 0 , 3845702E 06 EO V ECTOR = 0.37566 67E 00 

-0 . 11 26239E p9 -0.86 96515E 00 

-0 . 8864911E p8 -0.32Q2821E 00 


E VECTOR = 0 , 3756667E 00 V ' VECTOR = 0.3857833E 08 " OT£CTOR” 0T14T3X1 7*E~Dfe 

0.6884210E 00 0. 112670 5L 09 0.4658943E 05 

0.6204440E 00 “ ' ~ "b"7 a "699 356E Te 673444414E 06 


SIGMA 


0.4215347E 02 


B MATRIX BEFORE INTEGRATION 


0.52741E.O1 
0.65169E 01 
0.19S09E 02 
0.37237E o?. 
0.55233E 02 
0.65833 E 02 

0 7777oE 0 2 
0.774 0 7E 02 
0.68976E 02 
0.54487E 02 
0.37199F 02 
0.20 ? 55E 02 
0.62348E 01 
0.13754E 01 


0.29119E oo 
0.75789E 0 1~ 
0.21417E 0 2 
0.39077E 0 2 
0.56^10 e 02 
O.7O960E 02 
0.78119E 0 2 
0,769(j 7E (j2 
0.677&3E 02 

0.52824E 0 2 
0.35459E 0 2 
0.19285E U 2 

U.720U1E 01 
0.10157E ol 


0.64657E 00 
0 ;»7o68E 0 r 
0.23064E 0 2 
0.40918E o2 
0,tiB54PE 02 

0.7?032e 02 
0.78382E 0 2 
0.76326E 02 
0.6647J.E 02 

0.51137E 02 

0,3373()E 02 
0.1785BE 02 
0 . 63844 E ol 
0.71282E 00 


0.1099 0 E ol 
K58954E-oT" 

0.24746E o2 
0.42757E o2 
0.80144E 0 2 

D.73023e 02 
0 . 78559E 0 2 
-0.75S65E '02 
0.65134E 02 
0.49428E 02 

0 , 32ol8 £ 02 

0.16476E o2 
0.55482E ol 
0.4G557E '00 


0,1&4o 5E 
■BTlTPfBC- 
0.2&461E 
0 . 4459oE 
0.G1G87E 
0.73939 E 
0.78651E 
0.7492BE- 
0.63744E 
0.477 q2E 
0.30325E 
0.15142E 
0,47720^ 
0.27255E 


01 0 . 22639 E ol 

02 o71'246rE-u2- 

02 0.28211E 02 

oa o.RSRi3 E 02 

0 £ 0.63187E 02 

02 0.74779 E 02 

02 Q.78657E 0 2 

02-' o;74Tl-2E- 02 

02 0,&2305 £ 02 

02 0.45963E 02 

02 0.28G55E 02 

02 0.13857E 02 

Ol 0.4o56o E 01 
0 0' 0‘.'1322 R e 0 0 


0. 29669E p i 

0.?998l£ 02 
0.48222E 0 2' 
0.64G35E o£ 
0.75540E 02 
0.7B577E 0 2 
0.73223E- 02 

0,8o 820 £ 02 
0.44215E 02 
0,270ll E 02 
0.12623E 0 2 
0.34UU3E (ji 
0,42946 E_oi 


0.37448E ol 
"07T52TCE-02- 
U.31772E 02 
0,5ool3E 02 

y,GG02 &E 02 
0,7622l£ 02 
0.78HHE 02 
0.TZ263E 02 
U.59292E 02 
0,424f>0 £ 02 
0.25396E 02 
U.11443E 02 
0,28 U 48E ol 
O,20O29E.O2 


0.4b98yt pi U.55229E Ql 


0 .33b8iL 0 2 
0.51781E 0 2 
0.6/3&3E o2 

O.76820 t 02 
0.78l6o £ 02 
0. 7X233 1 '02' 
0 » 8 J72UE 0 2 
0,4o7u'4L 02 
0.23813E o2 
0.1U317E 02 
U,22fe9lt (ji 

o.otruffot. 'oo* 


0,354o4E 02 
0 . 5352 2 E q2 
U.68631E o2 

0,7733^ 02 
0.77825E 02 
"0,70 137 C-02' 
0.56122E o2 
0 , 3094-9E 02 
0.2226SE 02 

0.92476E ol 
0.17929E oi 

ro,X3366E-3B" 


B MATRIX AFTeR INTEGRATION 


O.OOOOOE 00 
0.259B7F 02 
0.1529 2F 03 
0.43&07E 03 
0 . 699&1 E 03 
0.15293E o4 
0.22738F 04 
0,3o568F 04 
0.3794°E 04 
0 . 4416 qE 04 
0.48752E 04 
O.51029E 04 
0.53036E 04 

'oTSWrer w 


U,1Gu91e 00 
u.33o3lE 0 2 
0.17364E 0 3 
0.47423E 0 3 

0.955&9E 0 3 

0.15997E 0 R 
0.23518E o4 

0.3134 0 E OR 
0.38632E ()4 
U.44&97E o 4 
0.49116E 04 
0.51829E 0 4 
0.53H3E o4 
0,53478r'o4 


0 , 62135fc, 00 
0.4H67E o2 
0.19577E 0 3 
0,t | l422 E 03 

0.10134E 0 R 
0,1S71PE a 4 
0.24300 e 04 
O.^PIOB^ o'R“ 

0.39304E 04 
0.R5217E 0 4 
0.494&2E 04 
0, s 20l5E o4 
0. 5 3lB2E o4 

0753>F87r-o4'" 


0.14664E 01 
0 . 5o465E 02 
U.21960E (j3 
0,556 q7E o3 
0.10727E 04 
0.17437E 04 
0,2Su85E q4 
0.328G6E oR 
0.39962E oR 
0.45719E o4 

0 , 4979qE o 4 
0.52l88 £ 04 
0.53241E q4 


THETO, CENTER Of INTENSITY ANGLE = 0.1777014E 02 

EPSILON. LOS ANGLE = 0.48B1530E 01 ARC-MINUTES 


0 ,28491 e 01 
0.6o98oE o2 
0.24528E £>3 
0.59974E 03 
0.11336E o4 
0.18172E 04 
0.25871E 04 
'07736l9E- 0 4 
0.40806E 04 
0 . 4&205E 04 
0.50102E 04 
0.52344E 04 
0.53293E o4 


0 , 47940£ 01 0 

0.72782E 02 (T 
0.?726l E 03 0 

0.64524E 03 0 

O.llSSlE 04 o 
0.18916E 04 0 

0.26658E 04 0 

ff734'3CTE” oR 0 

0.R123&E 04 0 

0.4G673E 04 0 

0.50397E 04 0 

0.52489E o4 0 
0.53337E o4 0 


,74036£ 01 
B592o E 0'2‘ 

,3{ll7lE o 3 
,69 2 56E o3 
,12GOO E U4 
.19667E oR 
.27444E oR 
.3?ToTE TPF 
.41052E oR 
,47l24f 04 
.50&75E o4 
.52622^ oR 
, 53374E y 4 


U.10753E 02 
0.TO0R5E 03 
U.332S8E q 3 
U.7416BE o3 
0.132&3E 04 
0 . 20 R 26 E oR 
U.28229E OR 
U . 35-829E-04' 
0.R2R53E OR 
U.47568E 04 
U , 5 q937E 04 
0.52742E 04 
0,b34obE y4 
~uT535'O0E 04’ 


O.1R910L 

0.11643E 

0.38b26E 

0.792588 

0.13920 E 

0.21192E 

0.29o12 £ 

' O'. 3 6 34 S t" 

0.4303 ut 
0,4 797it 
U,5ll«3t- 

0.528blt 

0.53430E 

"oToolTo 0*7 


02 0.19974l 02 
03 — 07T339XE o3- 

03 0 . 39975 E o 3 

03 0.84524E o3 

04 O.IRSOO^ 04 
OR 0 . 219 S 3 E 0 R 
OR 0 . 29792 E <)R 

■04 0.37233 E"o 4 — 

OR 0.436o7E OR 
04 0 . 48372 E o 4 

04 0 . 51414 E oR 

OR" 0 . 52949 E o 4 
04 0 . 53451 E 04 


EXECUflON TIME = 10 31,91 


- 

- 

- 

- 



- -- - 

: 





^ rv 









LISTING OF SUBROUTINE BIGNT 
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subroutine bignt, integral evaluator for lunar radiometric center 


c-eRRS...stno«c fortran source 


statements ioentfcn 


♦♦compiler messages^ 


CCCCCCCCCCLCCCCCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCC 

c 
c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

*c" 

c 


SUBROUTINE name - BIGNT 

utilization . dfTermines integration ranges 'aNd evaluates 
double integral for lunar Radiometric center function 


printed output 


UaTa SWITCH 13 ON PRINTS THETa-O, OELTa ALPHA 

and delta theta values. 

data SWITCH 14 ON PRINTS ALPHA VaLUES. 

DaTa SWITCH 0 ON PRINTS INTEGRATION ARRAYS BE- 
FORE and after first and second integrations. 

values input bv mainline - delta values for alpha and theta. 

DALP AND DTHET. 

VaLUls returned to mainline - aRRay of step-wise second integra- 
tion RESULTS. B > AND NUMBtR OF ELEMENTS IN ARRAY. U. 

ANGLES INPUT TO THIS SUBROUTINE ARE EXPECTED to be TO be IN DE- 
GREES. Values input to function 'flux* are in radians. 

SUBROUTINES REQUIRED - * FLUX * . INTEGRAND eVaLUaTO~R~ 

' QSF ' . SjMPSON * S RULE INTEGRATOR 
• LOGSW ’ . DATA SWITCH TESTER 


marda barthuli 
ball brothers research corporation 

BOULDER. COLORADO 


September 1974 


ll 




CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


SUBROUTINE BIGNT(DALP, DTHET, B, J) 
LOGICAL LOGSW 
COMMON SIGMA, RAD. UIAM 
DIMENSION 8(181), X( 90 ) 

DATA PI /3 ,141592654/, OELT/1 , 0E-06/ 


SINO(X) = SIN(X^PI/iao.O) 
COsD(X) = COS(X^PI/180.0) 
”TAN0T5fF = SlMOTxT/COSDfX ) ' 


CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C c 

c determine integration Range for theTa-o dependent on lunar c 

_C _ PHA s t angle. SIGMA! THEN SET U P MAIN THET a LOOP TEST VaLUE. J. C 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c C 

SIG = SIGMA^PI/160.0 . . 

IF (SIGMA - 90. U) 10,10,11 

10 THET = (SIGMA-90. 0) ♦ DTHET/2.0 


C 

11 


GO TO 12 

THET c OTHET/2.0 
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subroutine bignt, integral evaluator for lunar radiometric center 


C-ERRS. 4f STNO,C,,,.. FORTRAN SOUR C_ E S T A T t_ M_E_ N T S 




IDENTFCN 


♦♦COMPILER MEbSAGtS^ 


C 

12 


THETO s THET 

d = (90.0-THET)/DTHET +1 


DO 500 K=l,d 

IF(L0RSW(13>) WRITE(1,1000 > THETO 
THEO = THETO*PI/180.0 
DALPH = OALP 

C C 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c c 

C DETERMINE INTEbR A TlON RaNGL FOR aLPHa» DEPENDENT ON SIGMa AND C 

c t h eta-o. th^n set up inner auPh 6 loop test value, i. also c 

C ADJUST DELTA alpha TO FIT WITHIN THE RANGE EVLNLr, C 

c c 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


A = -TAND(SIGMA>*TAND(THET0) 
IF(ABSIA) - 1,0) 13,13,14 


l/) 

o 


c 

13 

c 

14 


c 

15 

c 

16 


c 

17 

C 

ia 


PI/2,0) 15,14,14 


A = 6-ACOS { A ) 
IF((A+DELT) ■ 

ALPH = 0.0 
Q = 90.0 
GO TO 18 


IF<SlND(SIGMA>*SIND(THETO) ) 16,17,17 
ALPH =0.0 

Q = A*180 . 0/PI 
GO TO 18 

ALPH = A*180 , 0/PI 
Q = 90,0 

ALPHA = ALPH 

I = <G-ALPHA+2*DELT> /DALPH 
DALPH = (Q-ALPHAI/I 
ALPH = ALPH + OALPH/2.0 
ALPHA = ALPH 

' ITTLWSWT 1 3 D VP I T'E 1 1 , 1 i" 6 o"l OA~ LTPFTJ dTHE f 


C C 

CCCCCCCCCCLCCCCCCCCCCrCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c c 

c evaluate integrand for the Range of alpha, then' integrate c’ 

C USING *QSF’« SIMPSON’S RULE SUBR OUTI NE, C_ 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c " “ c 

DO 100 L=1,I 

PHI = EASIN(COSD(THETO)*SIND(ALPHA) ) 

THETA = ATANlTAND(THETO)/COSDtALPHA) ) 


X(L > = FL0X1 PHI , TBETA TsIG.TFeO >“ 
IFIL0GSW(14)I WRITE! 1, 3000 ) ALPHA 
ALPHA = ALPH + DALPH»L 
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SUBROUTINE bIGNT* INTEGRAL EVALUATOR FOR LUNAR RAdIOHETRIc CENTER 


C-ERRS...STNO»C FORTRAN S- 0 U K C £ S A T L » E W T luE NlFcN »*(;OMPlLtH _ 


100 

C 


500 

C 


600 
C 


CONTINUE 

CALL QSR(DALPH*X,X,I) 

B ( K I = X(D*2.0 

THETO = THET + K*DTHET 

CONTINUE 

IM.NOT. LOGSW<0)> GO TO 650 
WRlTE(5,40OO) 

DO 6U0 KK=1,J,10 
MM = KK+9 

WRITE ( 5.4100) (B(LL), LL=KK,MM) 
CONTINUE 


CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCOCCCCCCCCCCCCCCCCCCCCC 


C 

c 

c 

c 

c 


PERFORM SECOND INTEGRATION ON UELTa THETa. BUNP THETa-0 and go 
DO IT AGAIN. RETURN TO mAInl 1 nE AFTER COMPLETIN' 3 ftAN^E OF 
THETA-0, 


ccccccccccLccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
c c 

650 

c 


CALL QSFIDTHETiB.B, U) 

IF ( .NOT . LOGSW(O)) RETURN 
WRITM 5*5000 ) 

DO 7U0 KK=1,J,10 
MM = KK+9 

WPITH 5.4100 ) (BILL) t LL=KK,MM) 

700 CONTINUE 
C 

RETURN 
C 
C 
C 
C 

1000 
1100 
2000 
3000 
4000 

"ijroD — FORi'ortr*"-*-, i or 
5000 FORMAT I * 0 * » * B MATRIX AFTER INTEGRATION*//) 
C 
c 
c 
c 


FORMAT ( ’THETO = *iE14.7) 

FOrMAtI ’DALPH = * »E14,7. • AND OtHET = *,E14.7) 
FORMAT I • I = *,I4, * AND J = *,I4) 

F0RMAf(5x ** alpha = *,ei4,t)' 

F0RMAT< •!’ , t B MATRIX BEFORE INTEGRATION*//) 


END 

VARIABLE ALLOCATIONS 
BLANK COMMON BLOCK 
SIGMA (R*6 c • =7FFD 


RA0(R*6 C)=7FFA 


DIAM(R*6 C)=7FF7 


EQUIVALENCES & INTERNAL VARIABLES 

_ pi 7 r*s — tswuv 


snf( n~*6 — rsTfircrs - 

THE0(R*6 )=000D OALPH ( )=0010 

S(R*6 )=001C ALPHA ( R*6 )=001F 


'THETTTR^S" VSWtTC - 

A ( K*6 ( =0013 

1 ( 1*2 ) = 0022 ' 


'THET OCR*? ) -001)9 

0EuT(R*6 >=0016 
L ( 1*2 7*0023 


KTT’FS TSTJWCT 

AlpH(k*6 )=0U19 

PHI(K*6 ) =0024 






LAUNCH AT KAUAI, HAWAII, 3 NOV. 1974, 1 OH 38M U.T. 

THET0 = -0.4734652E 02 

DAL PH = 0.1069554E 01 AND DTHET = 0.1000000E 01 

ALPHA = 0.5347773E 00 

ALPHA = 0 . 1 604332E 01 

ALPHA = 0 • 2673 886E 01 

ALPHA = 0 .374344 1 E 01 

ALPHA = 0*481 2996E 01 

ALPHA = 0 • 5 882550E 01 

ALPHA = 0.6952I05E 01 

THET0 = - 0 . 4634652E 02 

DALPH = 0*1 022502E 01 AND DTHET = 0.1000000E 01 

ALPHA. = 0*511 25 13E 00 

ALPHA = 0 • 1 533754E Oi 

ALPHA = 0 • 2556256E 01 

THET0 = -0 » 4534652E 02 

DALPH = 0. 1-026583E 01 AND DTHET « 0.1000000E 01 

ALPHA = 0*51 329 1 6E 00 

ALPHA = 0*1 539874E 01 

ALPHA = 0.2566458E 01 

ALPHA = 0 • 3593041 E 01 

THET0 = -0 » 4434652E 02 

DALPH = 0.1 0284 1 2E 01 AND DTHET = 0.1000000E 01 

ALPHA = 0.51420 60E 00 

ALPHA = 0.1 5426 18E 01 

ALPHA = 0.2571030E 01 

ALPHA = 0.21 08244E 02 

ALPHA = 0.221 I085E 02 

THET0 = - 0 . 4334652E 02 

DALPH a 0 . 1 009642E 01 AND DTHET = 0.1000000E 01 

ALPHA = 0. 50482 14E 00 

THET0 = -0 . 4234652E 02 

DALPH = 0.1011 884E 01 AND DTHET = 0.1000000E 01 

ALPHA = 0 . 50 59420E 00 

ALPHA « 0. 1517826E 01 

ALPHA = 0.2529710E 01 

ALPHA = 0 • 1 264855E 02 

ALPHA = 0 » 1 366043E 02 

ALPHA = 0.1 467231 E 02 

ALPHA = 0 . 2074362E 02 

ALPHA = 0.21 75550E 02 

ALPHA = 0 • 2276739E 02 

ALPHA = 0.3-187434E 02 

ALPHA = 0 . 3288623E 02 

EXAMPLE OF INTERMEDIATE TELETYPE OUTPUT 
FROM SUBROUTINE BIGNT 
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LISTINGS OF FUNCTION SUBPROGRAMS 
FLUX AND RDFCN 
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bbRc 


FUNCTION SUBPROGRAM FLUX, LUNAR BRIGHTNESS FUNCTION EVaLUaTOR 


C-ERRS...STNO*C FORTRAN SOUR CE ST A T t M E N T S 




ioentfcn **compilek MESSAGES** 


c ? 

ccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccclcccc 

c 

SUBROUTINE name - FLUX 


UflLiZATION «, EVaI-UaTES LUNaR BRIGHTNESS FUNCTION aT a' SINGLE 
POINT 

SUBROUTINES REQUIRED - 'RDFCN't LUNAR ALBEDO FUNCTION 
MAROA BARTHULI 

BALL BROTHERS RESEARCH CORPORATION 
BOULDER. COLORADO SEPTEMBER 1974 


c - 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C 

FUNCTION FLUX(PHI.THETA.SIG.THEO) 


c . 

C 

c 

C 

C 

C 

C 

C 

c 

c 


in 


COMMON SIGMA, KAD« 01AM 

CPHI = cos (PHI) 

CTHET = cos (THETA) 

CXI = CPHI*CTHET __ 

CGAM = CPHI*COS(THETA-SIG) 

Q = RAD/DIAM 

FLUX = (R0FCN<CGAM«CXI)*C0S<THE01 )/<1.0 + Q**2-2,0*CPHI*CTHET*Q) 
RETURN 


END 

VARIABLE ALLOCATIONS 
BLANK COMMON BLOCK 
SIGMa ( R*6 C)=TfFD 


RAD(R*6 C)=7FFA 


EQUIVALENCES s internal-variables 


FLUX ( R*6 
0 ( R*6 


) =0000 
>=OOOF 


CPHI(R*6 )=000 3 


DlAM(R*G C)=7FF7 
CTHET ( R*g ) =000& 


CXI(R*fc >=0009 


FEATURES SUPPORTED 
ONE WOR Q IN TEG ERS 
EXTENDED PRECISION 

CALLEO SUBPROGRAMS 
EcOS RdFcN EAUD 

REAL CONSTANTS 

.looooooooE or=ooiE 


Q 9 

^g 

I| 

*> > 

h—u M 


CGAM(K*6 )=oqoc 


INTEGER CONSTANTS 
2=0024 

CORE REQUIREM ENTS FOR 

BUAnk common 1 ' M» 


ESUB EMPY EdIV 
- ‘ ", 2000000 UOE 'o'l = du21 

FLUX 


ELd 


ESTO 


ESbR 


edvr 


EaXI SUbIN 


-ytn CONSTANTS AnD 1 PROGRAM- w 
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BBRC 


function subprogram rdfcn, lunar albedo function - Lambertian 


( 


C-ERRS.,,STNO.C FORTRAN S 0 U R C E_ S T A T t M E N_T S 


IOENTFCN 


**COMPILtH MESSAGES** 


C c 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCC 

c c 

c SUBROUTINE name - RDFCN > C 

c _ c 

c utilization - sets up lunaR albedo function, this function m a y c 
C Assume THE MOON to BE any kind oP radiator desired c 

C SUCH AS Lambertian, lommel-SelligeR, or combi- c 

c nations of the two. c 

C C 

C MARCIA BARTHULI C 

C ball BROTHERS RESEARCH CORPORATION c 


BOULOER. COLORADO 


September 1974 


cccccccccclccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 
C C 

FUNCTION RDFCN<CGAM,CXI) 
c 

COMMON SIGMA, RAD. DIAM 
C 

RDFCN = CGAMfcCXI 

return 


END 

VARIABLE allocations 
blank common block 

SIGMA(R*6 C>=TFFD 


rad < r*g c)=7ffa 


Ol AM ( R*6 C)=7FF7 


equivalences s internal variables 

RDFCN ( R*6 )=0000 

FEATURES SUPPORTED 
ONE WORD INTEGERS 
extenpfd precision 

"called 'subprograms 
empy ELD ESTO SUbIN 


CORE REQUIREMENTS FOR 

blank common- i°» 


rofcn 

variables and temporaries- 


constants a^d program- 


relative ENTRY POINT 1 ADDRESS lS“00O4 (HEX) 
END OF SUCCESSFUL COMPILATION 


COPYRIGHT UNA ^TemS , ' INC, 



LISTING OF 

INTERMEDIATE MAINLINE PROGRAM HERB 
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bbkc 


program herb* sigma inckementer for lunar center problem 


C-ERRS* FORTRAN 


SOURCE__SJ^TEMENTS IDENTFc_N_ _**COMPILEK MESSAGES** 


04 

oo 


C C 

CCCCCCCLCLLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCLCCCCC 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


PROGRAM NAME - HERB 

PROGRAM UTILIZATION 
PROBLEM, 


mainline' for luNaR Radiometric - center 


after reading in a delta sigma value <sigma c phase a ngle>, 
herb RUNS sigma FROM 0 TO 100 DEGREES* INTEGRATING THE 
BRIGHTNESS FUNCTION For EACH VAlUE AnD PLOTTING the re- 
sulting curve in terms of analytical angle psi and btosht- 

NESS. 

THE CENTER of BRIGHTNESS ANGLE IS CALCULATED FOR E A CH V A LUE OF 

sigma a nd printed along with the brightness Ra t i° of fulcto 
half moon. 

maroa barthuli 

ball BROTHERS RESEARCH CORPORATION 
BOULDER* COLORKaDO 10 SEPTEMBER 1974 


CCC CCCC CCC LCCC CCCC CCCCCC CCCCCCCCCCCLCCCCCC CCCC CCC CCCCCCCCCCCCCCCCCCCCC CC 

common sigma, rad* oiam 


OIMENSION Btieo). BB1 ( 19 > 
DAT* PI/3. 1H J59265 A/ 


"cccccccccclccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c _ c_ 

C INITIALIZATION section c 

C INITIALIZE LOOP values and set up plotting format c 

c ■ ' C 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

“c ~ - . - - — c 


Q © 

^ to 


SI 


c? 


RAO s 1736.0 
OIAM = 364403.0 

- REAO(2,1100) DALPH, DTHET, DSIG 


KK = 180 . OVb'sIG + T~ 

JJ = (180.0/DSlG)/2.0 + 1 


WR1TE( 1,4100 ) 

CALL SCALE< ,05*9.0,-90,0,0.0) 
CALL EGRIOIO, -90, 0,0. 0*10. 0,18) 
CALL EGRID(<S,0.0*l,0,.l,iO) 


C _ _ .C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c c 

C MAIN SIGMA LOOP ... DO INTEGRATIONS C 

c c 

ccccccccccCcccc'ccccccccccccccccccccccccccccccccccccccccccctcCcccCCLCcccc 


SIGMA = 0,0 
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bbRc 


program herb* sigma incRementeR for lunar center problem 


C-ERRS,.. STN P-»- C -**-t-*-L F 0 R_ T R A N S 0 U R C E _s 2_A__T N.J[ S__. 

DO 300 II-l , KK 

CALL BIGNT<DALPH,OTHET,B*U> 

BBK II ) = B(J) 

IF ( I I .EQ. <JJ> RATIO = ( BB1( II ) /BB1 ( 1) )*100 . 0 


IDENTFCN 


♦♦COMPILER MESSAGES** 


C „ 


CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c p 

c evaluate center of brightness angle for one Value of sigma c 

r C 


CCCCCCCLCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c - c 

HALF = 

THET = SIGMA - 90.0 

IF (SIGMA .GT. 90.0) THET = 0.0 

THETA = THET + DTHET/2, 0 


A 

l 

t>l 


00 200 1=1, J 

IF(B(I) - HALF) 101,202.203 
A1 = THETA 

THETA = THET + DTHET/2. 0 + I*OTHET 

.continue . . . 

AO = THETA 
GO TO 210 

A2 = THETA , „ „ . 

A <L = (HAL F *( A2- A1 )_ +_Bj I J|* A1 - B ( 1-1 )*A2 ) / ( B ( I ) - B(l-l)) 


C 

101 
C 

200 
C 

202 
C 

203 
C" 

C ... ...... --Sr 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c - - -- £ 

C PLOT ONE CURVE FOR ONE VALUE OF SIGMA C 


‘CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

C . . - - _ — c 

210 WRITEf 1,1)200) SIGMA, BBl(II). HALF, AO 
THET = SIGMA - 90.0 ___ 

IF (SIGMA .GTi 90.0) THET =0,0 

THETA = THET + DTHET/2. 0 


250 


DO 250 UJJ = 1,J 

CALL EPL0T(-2.THETA, B ( J JJ > /BB1 ( 1 ) > 
THETA = THET + DTHET/2. 0 + JJd*DTHET 
CONTINUE 

CALL UP 


C 

CCCCCCCCCCLCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c . . . c 

C GO do INTEGRATIONS FOR ANOTHER VAL0e OF SIGMA C 

C WHEN done WITH ALL THE SIGMAS EXCEPT SIGMA = ISO , WRITE T H E HALF C 

~C TO FULL'MOOKflJRIGHTNESS k PREtTT THEPlOT^ THEN GO' C 

C AWAY. c 

r • C 


£> 
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bbRc 


program herb* sigma incRementer for lunar center problem 


C"ERRS,,.stno_.c FOR t_r A_N_ SOU RCE STATEMENT S 


ioentfcn 


♦♦COMPILER MESSAGES** 


5200 FORMAT (F5,l) 

c 

5300 FORMAT(»BRIGHTNESS* > 

c __ 

5400 FORMAT ( F3 , 1 J 

c 

END 

variable allocations 
blank common block 

SISMa ( R*6 C ,=7 P FD RAD(R*G C>=7fFA 

EQUIVALENCES & INTERNAL variables 
DALPH(R* 6 >=0 2 00 OTHtT ( R*6 >=0 2 03 

11(1*2 J=020B J(I*2 )=020C 

THET A ( R*6 )=0216 I ( 1*2 >=3219 

J JJ ( 1*2 ) =0223 XX < R*G )=0224 

R(R*6 )=045A-0241 Bfil(R*6 (=0493-0450 


DIaM(R*6 C>=7fF7 


DSIG ( R*6 >=0 2 06 

rAtIo (K*6 )=0 20 0 

Ai(R*6 ) =021 A 
Y(R*6 >=0227 


KK ( 1*2 >=0 2 C9 

HAlF(R*G >=0210 
Ao ( R*S )=02lO 
X(R*6 )=022 a 


jju* 2 >=0aoA 
THE-j(R* 6 >=0213 
A2(R*6 )=02'2o‘ 

Pit K*6 >=022 d 


STATEMENT ALLOCATIONS 

1130=o4o9 4ioo=04dc 4200=0 <t FF 4300=o5o5 SooOrCSle 5100=0533 5200=0537 5300=0539 5400=0*40 

200=0629 202=0632 2qP=0 6 38 210=o66A 25o=O 60 E 300=q6D9 4oo=06E2 2o5=q73F 2u6=q77D 


10X-0&14 


FEATURES SUPPORTED 
ONE WORD INTEGERS 
EXTENDED PRECISION 
ORIGIN 
IOCS- 

PL OTTER 

1403 PRINTER" 

TYPEWRITER 

CARD 


CALLED SUBPROGRAMS 

Scale egrid siwr eplot up echaR init* eadd esub esubx empy 

EST’oX ESBR E'DVR WRTTZ "caRdZ"' VcRKT SR£0 SWR1 STOMP SFTd ST0F> 

IFIX FLOAT ISUb lrgt lreq 


EoIV tLD ELpX ESTO 

"SI OF sogstr “PRNZ SNR" 


real constants 

',173800000E 64=0496 
,900000000£ o1=04A5 


.100000000L 0 0“04B4 
,69(1000000 E 02= 04C3 


".3a4463"0OOE 06=5499' 
^OOOOOOOOE o2=q4A8 




,500000000E Ol=04C6 


.18'OOQOOOOE 03=049c 
OOOOOOOOOE 00=04AB 




.990000000E 00=0 4C9 


.'ffOOOOOOUOE 0T=049 f 
IOOOOOOOqE 02 = 04AE 




.700000000E 01=04CC 


"•500O0U06uE»01=04a 2 
.lOOOOUOOOE Ol=o4Bl 


Tfmmmiiimi-* 


,500O0U0(iljL U 1) =0 4CF 


INTEGER CONSTANTS 

2=0402 1=0403 


0=o4o4 


18=04D5 


3=04d6 


10=04d7 


7=04d8 


CORE requirements FOR 
Blank common- 10 * 


variables a«d temporaries- 662 » constants and program- 754 


end of SUCCESSFUL COMPILATION 


COPYRIGHT UNa SYSTEMS, INC 



EXAMPLE OF PLOTTED OUTPUT 
FROM MAINLINE PROGRAM HERB 
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Section V 

SYNOPSIS OF OPERATIONS 


The following steps gather, in logical order, the rather rambl- 
ing developments of the preceding sections. By knowing the 
location of the sensor (three space coordinates and the time) 
one can, by following the steps given below, compute the 
angular offset between the center of the moon and the lunar 
sensor's optical axis. 

These steps are, incidentally, precisely those performed by 
the computer program which is discussed in Section 4. 

1. Compute the right ascension and declination of the 
sun and moon using ephemeris data and Equation 4.2. 

2. Compute vectors G, H, and L (see ’Figure 3-9) using 
Equations 3.38, 3.41, and 3.39. The magnitudes of 

H and L are obtained from ephemeris data and the ap- 
plication of Equations 4,1 and 4.2. 

3. From G, H, and L, compute vector D using 3.40 and 
vector S* using 3.42. 


5-1 




4. Compute the phase angle, a, from D and S* using 3.37. 
If o<0, note that the center of brightness will be to 
the west of the selenographic center and use the ab- 
solute value of a in all subsequent calculations. 

5. Now compute limits of integration for 0o. 

5a. If o< tt/2, then a-ir/2 < 0 O < tt/2 
5b. If o>t r/2, then 0 < 6 0 < tt/2 

6. Now compute limits of integration for a. Evaluate 
3.26 to get do. The range of integration ove.r a 
is as follows: 


6a . 

If 

5a (above) applies, 

then 


• 

0 < a < a 0 

if 

O 

V 

o 

CD 


• 

0 < a < tt/2 

if 

0o > 0 

6b . 

If 

5b (above) applies , 

then 


• 

oto < a < tt/2 

if 

0o < cr-ir/2 


• 

0 < a < tt/ 2 

if 

0o > o-ir/2 
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7. Evaluate the integral appearing on page 3-23 over the 
ranges of 0o and a obtained in steps 5 and 6. The 
result is a table of values of this integral as 0 O 
spans its range. 

8. Using this table of values, determine (using 4.2) 
what value of 0o produces half the value produced by 

A 

00 = 11/2. Call this quantity 0 O . 

A 

9. Compute the sensor offset angle, e, using 0 O in 
Equation 3.44. If the original phase angle, a, 
(computed in step 4) was negative, negate e. 
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Appendix A 

Mathematical Notation 




Coordinate systems are designated by the letter E and are dis- 
tinguished from each other by superscripts. Thus, a system 
aligned to the celestial sphere could be called E - while one 
attached to the earth could be E e . Base vectors of coordinate 
systems are not termed x, y, and z but rather 1, 2, and 3 and 
form a right-handed set in that order. If confusion between 
systems is likely to occur, base vectors are referred to by 
subscripts on the designator of the system of which they are a 
part. Thus, the 2-axis of the earth-based coordinate system is 
termed E^. 

Vectors are invariably -represented as 3 x 1 matrices. Thus, a 
vector which would, in the notation of elementary vector analy- 
sis, be written as 


V = 3i - 4j + 7F 


is now expressed as 


V = 


3 

-4 

7 


A- 1 



The reason for this is that virtually all the vector operations 
that will be performed are scalar products and orthonormal ro- 
tations and the matrix formulation is just more convenient. 

The scalar, or "dot", product of two vectors (for example, A 

T 

and B) is written synonomously as A*B or A B where T denotes 
matrix transpose. Thus 


[ a l a Z a 3] 


T 

A*B = 


= Tj a j b j 

j=l 


L 3 -* 


Since sines and cosines are used so extensively in the defini- 
tion of vectors and rotation matrices, these functions have 
been abbreviated to: 


sine (Y) = sY 
cosine (Y) = cY 


A- 2 




Appendix B 

Something to Read if the Mathematics are Dull 



Another time, we went to Hannheim and 

ATTENDED A SHIVAREE - OTHERWISE AN OPERA - 
THE ONE CALLED "LOHENGRIN". THE BANGING AND 
SLAMMING AND BOOMING AND CRASHING WERE 
SOMETHING BEYOND BELIEF. THE RACKING AND 
PITILESS PAIN OF IT REMAINS STORED UP IN 
MY MEMORY ALONGSIDE THE MEMORY OF THE TIME 
THAT I HAD MY TEETH FIXED. "THERE WERE 
CIRCUMSTANCES WHICH MADE IT NECESSARY FOR 
ME TO STAY THROUGH THE FOUR HOURS TO THE 
END; AND I STAYED! BUT THE RECOLLECTION OF 
THAT LONG; DRAGGING; RELENTLESS SEASON OF 
SUFFERING IS INDESTRUCTIBLE. To HAVE TO 
ENDURE IT IN SILENCE; AND SITTING STILL; 

MADE IT ALL THE HARDER. I WAS IN A RAILED 
COMPARTMENT WITH EIGHT OR TEN STRANGERS; OF 
THE TWO SEXES; AND THIS COMPELLED REPRESSION; 
YET AT TIMES THE PAIN WAS SO EXQUISITE THAT 
I COULD HARDLY KEEP THE TEARS BACK. At THOSE 
TIMES; AS THE HOWLINGS AND WAILINGS AND 
SHRIEKINGS OF THE SINGERS; AND THE RAGINGS 
AND ROARINGS AND EXPLOSIONS OF THE VAST 
ORCHESTRA ROSE HIGHER AND HIGHER, AND WILDER 
AND WILDER, AND FIERCER AND FIERCER, I 
COULD HAVE CRIED IF I HAD BEEN ALONE, THOSE 
STRANGERS WOULD NOT HAVE BEEN SURPRISED 
TO SEE A MAN DO SUCH A THING WHO WAS BEING 
GRADUALLY SKINNED, BUT THEY WOULD HAVE 
MARVELLED AT IT HERE, AND MADE REMARKS ABOUT 
IT NO DOUBT, WHEREAS THERE WAS NOTHING IN 
THE PRESENT CASE WHICH WAS AN ADVANTAGE OVER 
BEING SKINNED. 


Each sang his indictive narrative in turn, 

ACCOMPANIED BY THE WHOLE ORCHESTRA OF SIXTY 
INSTRUMENTS; AND WHEN THIS HAD CONTINUED 
FOR SOME TIME, AND ONE WAS HOPING THEY MIGHT 
COME TO AN UNDERSTANDING AND MODIFY THE NOISE, 
A GREAT CHORUS COMPOSED ENTIRELY OF MANIACS 
WOULD SUDDENLY BREAK FORTH, AND THEN DURING 
TWO MINUTES, AND SOMETIMES THREE, I LIVED OVER 
AGAIN ALL THAT I HAD SUFFERED THE TIME THE 
ORPHAN ASYLUM BURNED DOWN. 

I HAVE SINCE FOUND OUT THAT THERE IS NOTHING 

the Germans like so much as an opera. They 

LIKE IT, NOT IN A MILD AND MODERATE WAY, BUT 
WITH THEIR WHOLE HEARTS. THIS IS A LEGITIMATE 
RESULT OF HABIT AND EDUCATION. OUR NATION 
WILL LIKE THE OPERA, TOO, BY-AND-BY, NO DOUBT. 

One in fifty of those who attend our opera 

LIKES IT ALREADY, PERHAPS, BUT I THINK A GOOD 
MANY OF THE OTHER FORTY-NINE GO IN ORDER TO 
LEARN TO LIKE IT, AND THE REST IN ORDER TO BE 
ABLE TO TALK KNOWINGLY ABOUT IT. . THE LATTER 
USUALLY HUM THE AIRS WHILE THEY ARE BEING 
SUNG, SO THAT THEIR NEIGHBORS MAY PERCEIVE 
THAT THEY HAVE BEEN TO OPERAS BEFORE. THE 
FUNERALS OF THESE DO NOT OCCUR OFTEN ENOUGH. 







Appendix C 

Evaluation o£ Lambertian and Lommel-Seeliger 
Reflectance Functions 




In this appendix, we wish to compute the half-moon to full- 
moon brightness ratios for three reflectance functions. 

The power density of the moon-reflected sunlight at some ob- 
servation point Is 


B = k J £( 'X^) dA 

A 


C.l 


where k 
R 
A 

f 

Y 

e 


a proportionality factor 

distance from the observation point to dA 
illuminated portion of the moon visible from the 
observation point 

lunar-surface reflectance function 
incidence angle of sunlight at dA 
reflectance angle of sunlight from dA to the ob- 
servation point 


If we make the assumption that the observation point to moon 
distance is very large compared to the lunar radius , we can 
treat R as constant and the visible region of the moon becomes 
a complete hemisphere. 


V 

V 


C-l 



Under these simplified condition, equation C.l can be written 
(in terms of normal spherical coordinates <f> and 0) as 


BCcr) 



a- tt/2 


it/ 2 

f £(y,5) s<j> d<j> de 
0 


C. 2 


where r = lunar radius 

a = lunar phase angle 

Our approach will be to evaluate C.2 for a = 0 (full moon) * 
and cr = tt/2 (half moon) with a particular function f, ratio 
the two terms and compare the result to the experimentally 
obtained value of 0.089 (ref. "Astronomical Quantities"; Sec- 
tion 3.4). 

Lambertain Reflector 

The reflectance function for a Lambertain reflector has the 
form 


f(y,£) = cy c £ 


C. 3 


C-2 



Now equations 3.8 and 3.19 give as 


C E , = ccj) C0 


C . 4 


CY = c<f> c( 8 - 0 ) 


C. 5 


and C.2 becomes 


tt/2 ir /2 


B(o) = 2k 




c <f> c9 c( 0 -a) s<f» d 4 > d 0 


a-ir/2 0 


C . 6 


which can be separated into 


tt/2 


tt/2 


B(a) = 2 k { ^ ) / c 1 s<j> d<{> J c 0 c( 0 -o) d 0 
0 a-ir /2 


C . 7 


Since the <j> - integral is simply a constant, it can be absorbed, 

2 

together with the 2k(r/R) term into a single super-constant, 


C-3 




K, which leaves 


tt/2 

B (a) = K fee c(0-a) de 


C . 8 


a-Tr/2 


This can be integrated directly to give 


B 


Ca) = — |(tt-o)ccf + so | 


C.9 


Thus 


B (tt/2) 
“B CO) 


= 0.32 


C.1C 


Lommel-Seeliger Reflector 

The reflectance function for a Lommel-Seeliger reflector has 
the form 


£ (y »0 


cy 

cy + c£ 


C. 1] 


C-4 



which, by using C.4 and C.5, transform into 


f(Y,S) 


c(9-a) 
c ( 9- a ) + c0 


C.l 


and 


BO) 



c(tOOe d * « 


a~ ir/2 0 


C.l 


As before, the <j> - integral can be extracted- and constant terms 
combined. Doing so yields 

ir/2 

" K f c(| - -ar l C6 de 

c-tt/2 

Rather than attempt to integrate C.13 as it stands, let us 
instead perform the integration twice; once with a = 0 and 
once with a = ir/2. 


C-5 




Thus , 


B COD = K / = Ktt/2 


C. 15 


-ir/2 


5 • 


- In | c8 - s0 


= Ktt/4 


C.li 


B C TT/ 2 ^ _ j-v r 

— Trrm — u * D 


C-6 



Something Else 


In the discussions of Section 3.4, it was stated that a "modi- 
fied Lommel-Seeliger" reflectance function was also investigated. 
It had the form 


fCY.O 


c 2 Y 


cy + c£ 


C. 20 


Putting this in terms of <f> and 6 gives 


f CY»0 


c<j>c 2 (8-g) 
c(e-o] + c9 


C. 21 


and the resultant brightness function , B(cr) , becomes 


B(c) 


it/ 2 



c 2 (8-q) 
c(0-cr) + ce 


d0 


o-ir/2 


C . 22 


C-7 



As, before, we compute B(0) and B(tt/ 2) separately: 


B(0) = 


ir/2 

j I c6 d9 


C . 2 3 


it/ 2 


= K 


C , 24 


ir/2 


B CW2) - *l*rrcs ds 


C. 25 


0 


The author spent his normal, self-imposed time allotment of 5 
minutes trying to integrate C.25 before turning to his best 
beloved table of integrals. After another 10 minute search, 
he abandoned that approach as well in favor of a numerical 
solution using Simpson’s rule. The integrand was evaluated at 
eleven points and an error of less than 10 ^ was expected. 


C-8 



The result of this numerical integration gave 


B(ir/2) = 0.62 k 


C . 26 


so 


B 0/2) 

B ( 0 ) 


0.62 


C . 2 7 


C-9 




Appendix D 

Predicted Sensor Offsets 
for 

Kauai and White Sands Launches 


Taken from the letter: 

"Radiometric Center of the Moon Study" 
Charlie Rose to Morris Gisser 
20 September 1974 




The sensor offset angle, e, has been calculated for two launches; 
Kauai, Hawaii and White Sands, New Mexico. The offsets were 
calculated for the opening, center, and closing of the launch 
window in each case. 


Location 

Time 

CU.TO 

Offset Angle, e 

Kauai, Hawaii 

10 h 

38 m 

4.88 

arc minutes 

3 November 1974 

ll h 

20™ 

4.92 

If If 


12 h 

3 m 

4.98 

It ?! 

WSMR 

5 h 

3 m 

1.5-5 

arc minutes 

28 December 1974 

5 h 

41 m 

1.51 

it it 


6 h 

19 m 

1.47 

1! ft 


D-l 



